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PREFACE 


This  report  is  an  interim  product  of  a joint  effort  of  the  Department  of 
Health  and  Human  Services  and  the  National  Bureau  of  Standards  (NBS) , Center  for 
Fire  Research.  The  program  is  a multi-year  activity  initiated  in  1975.  It  con- 
sists of  projects  in  the  areas  of:  decision  analysis,  fire  and  smoke  detection, 

smoke  movement  and  control,  automatic  extinguishment,  and  behavior  of  institutional 
populations  in  fire  situations. 

This  report  describes  a computer  program  which  analyzes  pressurized  stairwells 
and  pressurized  elevators.  The  program  was  initially  intended  as  a research  tool 
to  investigate  the  feasibility  of  specific  systems.  However,  this  program  may  be 
of  interest  to  design  engineers  responsible  for  pressurized  stairwells  or  pressurized 
elevators . 
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Figure  1.  Flow  chart  for  main  program 


A COMPUTER  PROGRAM  FOR  ANALYSIS  OF  SMOKE  CONTROL  SYSTEMS 


John  H.  Klote 


Abstract 

This  paper  describes  a computer  program  developed  to 
analyze  systems  intended  to  control  smoke  in  building 
fires.  These  systems  include  pressurized  stairwells, 
pressurized  elevator  shafts,  zone  smoke  control  systems, 
and  pressurized  corridors.  This  program  calculates  air 
flows  and  differential  pressures  throughout  a building 
in  which  a smoke  control  system  is  operating.  The  basic 
assumptions  and  limitations  of  the  program  are  also  dis- 
cussed. The  appendices  contain  a program  listing  and 
examples . 

Key  words:  Air  movement;  computer  programs;  egress; 

elevator  shafts;  escape  means;  modeling;  pressurization; 
simulation;  smoke  control;  stairwells. 


1.  INTRODUCTION 

The  majority  of  fire  fatalities  result  from  smoke  inhalation.  As  a result  of 
this,  a number  of  systems  have  been  designed  and  built  to  control  smoke  movement  in 
building  fires.  The  most  common  smoke  control  systems  are  pressurized  stairwells 
and  zone  smoke  control  systems1 2.  These  systems  are  intended  to  control  smoke  move- 
ment in  a building  by  use  of  air  flows  and  by  differential  pressures.  The  computer 
program  described  in  this  paper  provides  a means  to  calculate  the  air  flows  and 
differential  pressures  throughout  a building  (either  real  or  conceptual)  in  which  a 
smoke  control  system  is  operating. 

A number  of  computer  programs  have  been  developed  which  are  applicable  to  smoke 

control.  Some  of  these  programs  calculate  steady  state  air  flow  and  pressures 

2 

throughout  a building  [1,2]  . Other  programs  go  beyond  this  to  calculate  smoke 
concentrations  throughout  a building  that  would  be  produced  in  the  event  of  a fire 
[3-7].  In  general,  most  of  these  programs  are  capable  of  analyzing  smoke  control 


1The  concept  of  extending  the  use  of  smoke  control  to  protect  elevators  is  currently 
being  investigated  at  NBS. 


2 

Numbers  in  brackets  refer  to  the  literature  references  listed  at  the  end  of  this 
paper. 
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systems.  However,  the  program  described  in  this  paper  has  been  specifically  written 
for  analysis  of  smoke  control  systems  and  is  an  extension  of  a program  specifically 
written  for  analysis  of  pressurized  stairwells  and  elevators  [8] . While  the  basic 
theory  of  this  program  is  the  same  as  that  of  the  stairwell  program  it  has  been 
extended  to  include  analysis  of  (1)  stairwells  with  vestibules,  (2)  elevators  with 
elevator  lobbies,  (3)  zone  smoke  control  systems,  and  (4)  pressurized  corridors. 

The  data  input  has  been  designed  to  minimize  the  quantity  of  required  data  and  still 
maintain  a high  level  of  generality  in  the  model.  The  output  consists  of  the 
pressure  differences  across  all  of  the  building  shafts,  as  well  as  the  flows  and 
pressures  throughout  the  building. 

This  program  was  originally  intended  primarily  as  a research  tool  to  investigate 
the  feasibility  of  specific  smoke  control  systems  and  to  determine  the  interaction 
between  these  systems  and  the  rest  of  the  building.  The  predecessor  [8]  of  this 
program  has  already  been  used  to  analyze  pressurized  stairwells  without  vestibules 
and  to  evaluate  factors  which  affect  the  performance  of  these  systems  [9] . And, 
this  program  has  been  used  to  generate  data  for  an  National  Bureau  of  Standards 
(NBS)  Handbook  on  Smoke  Control  Design  which  is  being  developed.  This  paper  is  not 
intended  to  be  a design  guide  for  smoke  control  systems.  The  state-of-the-art  of 
these  systems  is  still  under  development  and  designers  of  these  systems  should  seek 
the  most  current  data  available. 


2 . PROGRAM  CONCEPT 

In  this  computer  program  a building  is  represented  by  a network  of  spaces  or 
nodes  each  at  a specific  pressure  and  temperature.  The  stairwells  and  other  shafts 
are  modeled  by  a vertical  series  of  spaces,  one  for  each  floor.  Air  flows  through 
leakage  paths  from  regions  of  high  pressure  to  regions  of  low  pressure.  These 
leakage  paths  are  doors  and  windows  which  may  be  opened  or  closed.  Leakage  can  also 
occur  through  partitions,  floors,  exterior  walls  and  roofs.  The  air  flow  through  a 
leakage  path  is  a function  of  the  pressure  difference  across  the  leakage  path. 

In  this  computer  model  air  from  outside  the  building  can  be  introduced  by  a 
pressurization  system  into  any  level  of  a shaft  or  even  into  other  building  spaces. 

This  allows  simulation  of  stairwell  pressurization,  elevator  shaft  pressurization, 
stairwell  vestibule  pressurization,  and  pressurization  of  any  other  building  space. 

In  addition,  any  building  space  can  be  exhausted.  This  allows  analysis  of  zone 
smoke  control  systems  where  the  fire  zone  is  exhausted  and  other  zones  are  pressurized. 
The  pressures  throughout  the  building  and  flow  rates  through  all  the  flow  paths  are 
obtained  by  solving  the  air  flow  network  including  the  driving  forces  such  as  the 
wind,  the  pressurization  system  or  an  inside  to  outside  temperature  difference. 
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3.  ASSUMPTIONS  AND  LIMITATIONS 


1.  Each  space  is  considered  to  be  at  one  specific  pressure  and  one  specific 
temperature . 

2.  The  flows  and  leakage  paths  are  assumed  to  occur  at  mid-height  of  each  level. 

3.  The  net  air  supplied  by  the  air  handling  system  or  by  the  pressurization 

system  is  assumed  to  be  constant  and  independent  of  building  pressure. 

4.  The  outside  air  temperature  is  assumed  to  be  constant. 

5.  The  barometer  pressure  at  ground  level  is  assumed  to  be  standard  atmospheric 

3 

pressure  (101325  P ) . 
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4 .  EQUATIONS 


A.  Flow  equation 


F = CA  /2pAP 


(3.1) 


where : 

F = mass  flow  rate 
C = flow  coefficient 
A = flow  area 

p = density  of  air  in  flow  path 
AP  = pressure  difference  across  flow  path 


The  flow  coefficient  is  dimensionless  and  for  smoke  control  analysis  it  is 
generally  taken  to  be  in  the  range  of  0.6  to  0.7.  Because  of  the  large  number  of 
flow  calculations  performed  during  the  computer  analysis  the  flow  equation  is 
rewritten  in  the  program  as  F = C'  /AP.  Using  the  ideal  gas  law,  the  adjusted  flow 
coefficient,  C',  can  be  expressed  as 


C' 


CA  J 2 Patm~ 
1 RT 


where: 

Patm  = abs°lute  barometric  pressure  at  ground  level 
R = gas  constant  of  air 

T = absolute  temperature  of  air  in  flow  path 


(3.2) 


3 

The  results  of  the  program  are  not  very  sensitive  to  changes  in  atmospheric 
pressure.  For  altitudes  considerably  different  from  sea  level  the  more  accurate 
value  can  be  substituted  by  changing  an  assign  statement  in  the  subroutine  INPUT 
and  one  in  the  subroutine  CORR. 
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B . Mass  Balance  Equations 


For  building  compartment  i 


S F(i,j)  + 2 

j=l  k=l 


F + F 

o(i,k)  *f(i) 


0 


(3.3) 


and  for  shafts 


where : 


[i,  j)  + Fo(i, 


k)  + Ff(i) 


k=l 


0 


(3.4) 


F(i  i)  = mass  flow  rate  from  space  j to  space  i.  For  building 

' -1  compartments  this  flow  can  be  either  horizontal  or  vertical, 
however  for  shafts  this  flow  can  only  be  horizontal. 


F . i k)  = mass  flow  rate  from  direction  k outside  of  the  building 
' to  space  i. 

Ff,.,  = net  mass  flow  rate  of  air  due  to  the  air  handling  system 
1 or  due  to  a pressurization  system. 

Nc  = number  of  building  spaces  connected  to  space  i. 

Nq  = number  of  connections  to  the  outside  from  space  i. 

is  the  space  number  at  the  bottom  level  of  the  shaft  and  the  spaces  in  the 
shaft  are  numbered  consecutively  up  to  ^ which  is  the  space  number  at  the  top  of 
the  shaft. 


C.  Shaft  Pressures 


The  following  relationship  is  used  to  calculate  the  gage  pressure,  P^,  at 
floor  i of  a shaft  in  terms  of  P^_^  at  floor  i - 1. 


P...  = P , . . > -P 

(i)  (l-l)  z 


P 


f 


where : 

Pz  = hydrostatic  pressure  difference 
P^  = pressure  loss  due  to  friction 


(3.5) 


The  following  equation  is  used  to  calculate  the  hydrostatic  pressure  difference. 


4 

In  this  paper  the  term  building  compartment  refers  to  a space  in  a building  other 
than  in  a shaft. 
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p 

RT 


h(i)  h ( i-1) 


3.6 


where : 

= height  of  point  i 
h^_jj  = height  of  point  i-1 

g = gravitational  acceleration 
R = gas  constant 


T = 


T(i)  + T (i-1) 


P , . . + P , . , . 

P = (l)-  X ^ ---  + P. 

2.  D 


P^  is  a constant  used  to  convert  an  average  gage  pressure  to  the  average 
absolute  pressure,  P. 

The  following  equation  is  used  to  calculate  the  pressure  loss  due  to  friction, 

f.  \ 2 
Im  a 

Pf  = slcHl  (3.7) 


where : 


= upward  flow  from  i-1  to  i in  shaft 

C = shaft  flow  coefficient 
s 

S = sign  of  mu 
D . Outside  Pressures 

Outside  pressures  can  either  be  entered  by  the  user  or  can  be  calculated  by  the 
following  method. 


o (i) 


P,  , . . + C P , . . 
h ( l)  w v ( l) 


(3.8) 


where : 
P 


o (i) 


h ( i) 

'v(i) 

C 

w 


outside  gage  pressure  at  height  h(i)  above  absolute  pressure 
at  ground  level 

hydrostatic  pressure  difference  between  h(i)  and  ground  level 
velocity  pressure  due  to  the  wind  at  height  h(i) 
pressure  coefficient 


Because  the  outside  temperature  is  constant 


Ph(i)  * patm  exp 


I - qhU)  \ 
l RTout  ) 


- p. 


(3.9) 
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where : 

Patm  = absolute  barometric  pressure  at  ground  level 
TQut  = outside  absolute  temperature 

When  the  outside  pressures  are  calculated  by  the  computer  the  wind  velocities 
are  assumed  to  be  described  by  the  power  law. 


where: 

VQ  = wind  velocity  at  height  hQ 
n = wind  exponent 

This  relationship  has  been  extensively  used  to  describe  the  boundary-layer 
velocity  profile  of  the  wind  near  the  surface  of  the  earth.  It  assumes  that  the 
terrain  surrounding  the  building  is  homogeneous.  That  is,  that  there  are  no  large 
obstructions  near  the  building  which  could  produce  local  wind  effects.  A value  of 
0.16  for  the  wind  exponent  is  appropriate  for  flat  terrain.  The  wind  exponent 
increases  with  rougher  terrain.  For  very  rough  terrain  such  as  urban  areas  a value 
of  0.40  would  be  appropriate. 

The  equation  for  the  velocity  pressure  at  height  h(i)  is  obtained  by  substituting 
the  velocity  from  the  power  law  into  the  usual  relation  for  velocity  pressure  (P  = 

|pvM. 


where  p is  the  outside  air  density. 

The  pressure  coefficients  are  in  the  range  of  -0.8  to  0.8  where  positive 

values  are  for  windward  walls  and  negative  values  are  for  leeward  walls.  The  z 

pressure  coefficient  depends  upon  building  geometry  and  varies  locally  over  the  wall 

surface.  Numerical  values  for  C and  n as  well  as  practical  engineering  information 

w 

are  available  from  a number  of  sources  [10-13]. 

5.  PROGRAM  DESCRIPTION 

This  program  is  written  in  ANSI  FORTRAN  on  the  UNIVAC  1100/82  and  a program 
listing  is  provided  in  appendix  D.  The  following  is  a detailed  description  of  the 
main  program  and  the  major  subroutines. 
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5.1  Main  Program 

The  main  program  calls  the  subroutines  which  read  the  data,  calculate  the 
adjusted  flow  coefficients,  calculates  the  initial  values  of  pressures  and  inter- 
atively  solves  for  the  pressures  according  to  the  logic  illustrated  in  the  flow 
chart  of  figure  1. 

5.2  INPUT  Subroutine 

This  routine  reads  the  data  that  are  necessary  for  a flow  analysis  of  the 
stairwell  or  elevator,  including  an  analysis  of  the  rest  of  the  building.  These 
data  consist  of  the  following: 

1.  Outside  temperature. 

2.  Temperature  throughout  the  building. 

3.  Outside  pressures.  These  can  be  entered  or  calculated  as  described 
earlier . 

4.  Description  of  the  flow  network  including  flow  coefficients  and  flow  areas 
for  all  connections  and  the  net  air  flows  to  each  space  due  to  the  air 
conditioning  system  or  due  to  a pressurization  system. 

The  data  above  can  be  entered  in  either  SI  units  or  in  engineering  units. 
Appendix  A contains  a detailed  description  of  the  data  input  method. 

In  addition  to  reading  data,  this  subroutine  provides  temperature  and  pressure 
data  as  well  as  a complete  description  of  the  flow  network.  This  routine  also 
calculates  initial  estimates  of  the  hydrostatic  pressure  differences.  When  data  is 
entered  in  engineering  units  the  subroutine  UNITS  is  called  which  converts  all  units 
to  the  SI  system. 


5 . 3  CORR  Subroutine 

This  routine  calculates  adjusted  flow  coefficients  for  all  flow  paths  using  eq. 
(3.2).  Two  sets  of  these  coefficients  are  calculated  for  each  flow  path  to  allow 
for  flow  in  either  direction. 


5.4  INIT  Subroutine 

This  routine  calculates  initial  estimates  of  the  building  pressures  by  a technique 
used  by  Sander  [1] . In  this  technique,  mass  flows  are  considered  linear  functions 
of  differential  pressure  and  therefore  the  flow  equations  can  be  expressed  and 
solved  in  matrix  form.  In  this  estimate,  shaft  pressures  are  considered  hydrostatic. 
The  resulting  pressures  form  a starting  point  for  the  iterative  solution  which 
follows . 
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Figure  1.  Flow  chart  for  main  program 
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5.5  BLDGP  Subroutine 


The  iterative  solution  for  the  building  pressures  and  flows  consists  of  t e 
three  subroutines  BLDGP,  SHAFTP  and  PZAD.  The  subroutine  BLDGP  operates  on  the 
building  compartments  sequentially.  The  sum  of  all  the  mass  flows  into  compartment 
i is  calculated.  If  the  absolute  value  of  this  sum  is  less  than  a convergence  limit 
then  eq.  (3.3)  is  considered  satisfied  and  the  computer  proceeds  to  the  next  com- 
partment or  returns  to  the  main  program.  However,  if  the  absolute  value  of  the  sum 
is  greater  than  the  convergence  limit,  then  an  improved  estimate  of  the  pressure  at 
compartment  i is  obtained  by  the  regula  falsi  method  [14] . When  none  of  the  pres- 
sures need  to  be  modified  this  routine  passes  a convergence  signal  to  the  main 
program. 

5.6  SHAFTP  Subroutine 

The  structure  of  this  routine  is  very  similar  to  that  of  BLDGP  except  that  it 
operates  on  shafts  sequentially.  The  sum  of  all  the  mass  flows  into  shaft  i is 
calculated.  If  the  absolute  value  of  this  sum  is  less  than  the  convergence  limit 
then  eq.  (3.4)  is  also  considered  satisfied  and  the  computer  proceeds  to  the  next 
shaft  or  returns  to  the  main  program.  However,  if  the  absolute  value  of  the  sum  is 
greater  than  the  convergence  limit,  then  improved  estimates  of  the  shaft  pressure 
are  calculated.  This  is  done  by  changing  the  pressures  at  the  bottom  of  the  shaft 
and  then  recalculating  the  shaft  pressure  by  eq.  (3.5).  Again  the  regula  falsi 
method  is  used,  and  if  none  of  the  shaft  pressures  require  modification  a convergence 
signal  is  passed  to  the  main  program.  It  can  be  seen  from  figure  1 that  if  conver- 
gence is  achieved  in  both  BLDGP  and  SHAFTP,  then  the  subroutine  OUT  will  print  the 
solution.  Otherwise,  the  hydrostatic  pressure  differences  are  adjusted  in  the 
subroutine  PZAD. 

5.7  PZAD  Subroutine 

This  routine  calculates  hydrostatic  pressure  differences  by  eq.  (3.6)  using  the 
most  recent  pressure  estimates. 

5.8  OUT  Subroutine 

This  routine  outputs  mass  flows  and  pressures  for  the  flow  network  as  well  as 
the  differential  pressures  across  each  shaft.  If  the  data  input  was  in  engineering 
units  then  appropriate  variables  are  converted  to  the  engineering  system  before 
output. 

6.  FUTURE  DIRECTION 

It  is  planned  to  use  this  computer  program  in  a project  at  NBS  to  study  the 
feasibility  of  protected  elevators  as  a means  of  fire  escape  for  handicapped 
individuals.  Consideration  is  being  given  to  further  development  of  the  program  for 
use  as  a design  tool.  Also,  a program  may  be  developed  for  microcomputers  which  can 
be  used  interactively. 
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APPENDIX  A.  DATA  INPUT  DESCRIPTION  FOR 
COMPUTER  PROGRAM 

Data  input  consists  of  the  following  elements: 

1.  Initial  data 

2.  Building  heights 

3.  Temperature  profiles 

4.  Outside  pressure  profiles 

5.  Building  data 

6.  Shaft  data 

Each  of  these  input  elements  is  described  in  detail  in  the  following  sections. 
Elements  1 through  6 are  always  required.  In  the  following  sections  the  input 
required  for  each  of  the  six  data  elements  is  described  in  detail.  Each  block  or 
group  of  blocks  below  represent  an  input  card.  Unless  otherwise  stated  these  cards 
are  unformatted;  that  is,  the  numbers  do  not  have  to  be  placed  in  specific  columns 
and  integers  can  be  written  with  or  without  decimal  points.  However,  separate 
pieces  of  numerical  data  must  be  separated  by  one  or  more  spaces.  Examples  of 
input  data  are  provided  in  Appendix  B. 

1.  Initial  data 


project  title  (col.  1-72) 


outside  unit  indication  summary  output  ^ 

temperature  (°C,°F)  (1  for  SI,  2 for  Eng)  (0  for  none,  or  file  number) 1 


1 


The  user  must  assign  this  file  before  program  execution. 
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2 . Building  heights 


N,  , no.  of  building 
levels 


input  parameter 
(either  0 or  1) 


If  input  parameter  = 0,  then  heights  for  each  building  level  are  to  be  individually 
entered  as  follows: 


h(l)  h ( 2 ) 

□ □ 


where  h 


(i) 


is  the  height  of  the  center  of  level  i above  the  ground  (m,  ft) . 


If  input  parameter  = 1,  then  the  following  card  must  be  entered. 

h,.,  distance  between 

floors  (m,  ft) 


3 . Temperature  profiles 


no.  of  temperature  profiles 


For  each  temperature  profile  the  following  data  must  be  supplied. 


no.  of  temp.  level  temperature  level  temperature 
points  no.  (°C,°F)  no.  (0C,°F) 


level  temperature 
no.  (°C,°F) 
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4 . Outside  pressure  profiles 


no.  of  outside 
pressure  profiles 


input  parameter 
(either  0 or  1) 


If  the  input  parameter  = 0,  each  outside  pressure  profile  is  entered  as  follows: 


Po(l)  Po  (2)  Po(3)  Po(i)  P°(V 


where  P , . . 

o (l) 


is  the  outside  pressure  at  the  center  of  level  i. 


If  the  input  parameter  = 1,  the  outside  pressures  are  calculated  and  the  following 
data  are  required. 


Vo 

wind  velocity 
(mph) 


h 


o 

height  at 
velocity  is 


which 

measured 


n 

wind 

exponent 


pressure  coefficients  for  each  pressure  profile 


CW(1) 


CW(2) 


CW(Npo) 


o o o o 


5 . Building  data 


N 


f 


no.  of  levels 
(or  floors) 
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All  of  the  following  data  in  this  input  element  are  supplied  for  each  level,  or 
consecutive  groups  of  similar  levels. 


I 


1 


Starting 


floor 


I 


2 


Ending  floor 


N 

com 

No.  of  compartments  per  floor 


(Floor  data  is  entered  in  ascending  order  of  levels  or  floors.  When  data 
is  for  only  one  level  then  1^  = and  the  same  number  is  supplied  for 

both. ) 


For  each  compartment  on  a level  the  following  data  are  supplied: 


NCS 

No.  of  connections 
to  other  compart- 
ments on  the  same 
level 


NCA 

No.  of  connections 
to  compartments  on 
the  level  above 


NC0 

No.  of  connections 
to  the  outside 


Net  flow^ 
(1/s,  cfm) 


Temperature 

profile 

number 


For  each  connection  between  this  compartment  and  another  on  the  same  floor  the 
following  data  are  required. 


Other  compartment 
number  on  the  same 
level 


C 

flow  coefficient 


flow 
(m2 , 


A 


area 

ft2) 


For  each  connection  between  this  compartment  and  one  on  the  level  above  the  following 
data  are  required. 


2 


All  net  flows  are  at  standard  conditions  of  21°C  (70°F)  and  one  atmosphere. 
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Other  compartment 
number  on  floor 
above 


C A 

flow  coefficient  flow  area 

(m2,  ft2) 


For  each  connection  to  the  outside  the  following  data  are  required. 


outside  pressure 
profile  number 


C 

flow  coefficient 


A 


flow  area 
(m2,  ft2) 


6 . Shaft  data 

no.  of  shafts 


All  of  the  following  data  in  this  input  element  are  required  for  each  shaft, 
shaft  title  (col  1-20) 


shaft  flow 
coefficient 


bottom 

level  of  shaft 


top  level 
of  shaft 


temperature 
profile  number 


Enter  the  following  typical  data  which  applies  to  each  level  of  the  shaft. 


Exceptions  can  be  entered  later. 


15 


no.  of  connections 
between  typical 
level  of  shaft  and 
outside 


net  flow  into 
typical  level 
of  shaft 
(1/s,  cfm) 


The  connection  data  to  the  building  for  a typical  level  are  required. 


compartment  no. 
to  which  shaft  is 
connected 


C A 

flow  coefficient  flow  area 

(m2,  ft2) 


For  each  connection  to  the  outside,  the  connection  data  for  a typical  floor  are 

required . 


outside  pressure 
profile 


C 

flow  coefficient 


A 

flow 
(m2 , 


area 

ft2) 


The  number  of  exceptions  to  the  typical  data  is  required. 


no.  of  exceptions 


All  of  the  following  data  in  this  input  element  are  required  for  each  exception. 


exception  type  level  of  shaft 

(1,  2 or  3) 
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The  next  card  depends  on  the  exception  type.  For  exception  type  = ' . 
to  the  net  flow  into  the  floor  of  the  shaft  is  defined. 


net  flow 
(1/s,  cfm) 


For  exception  type  = 2 , an  exception  to  an  outside  connection  for  this  shaft  is 
defined . 


outside  pressure 
profile  number 


C 

flow  coefficient 


A 


flow  area 
(m2,  ft2) 


For  exception  type  = 3,  an  exception  to  the  connection  between  the  shaft  and  the 
building  is  defined. 


compartment  no. 
to  which  shaft 
is  connected 


C 

flow  coefficient 


A 

flow  area 
(m2,  ft2) 
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APPENDIX  B.  INPUT  EXAMPLES 


1.  Example  1 

A ten  story  building  with  a pressurized  stairwell  and  no  vertical  leakage 
within  the  building  is  heated  to  70°F  when  the  outside  temperature  is  -20°F.  The 
stairwell  temperature  is  60°F  at  the  tenth  floor  and  50°F  at  the  bottom  floor.  The 
stairwell  is  pressurized  by  a net  550  cfm1  per  floor.  The  wind  is  30  mph  at  a 
height  of  30  ft  and  the  wind  exponent  is  0.14.  This  building  has  connections  to 
the  outside  in  two  directions.  The  wind  pressure  coefficients  are  0.7  for  the 
windward  wall  and  -0.7  for  the  leaward  wall.  The  flow  areas  are  the  same  vertically 
and  are  listed  in  Table  Bl.  The  flow  coefficient  is  taken  to  be  0.65  for  all 
connections . 


Table  Bl.  Flow  areas  for  example  1 

Connection  location  Area  (ft2) 

Between  stairwell  & building  0.42 

Between  building  & outside  into 

the  wind  0.75 

Between  building  & outside  away 

from  the  wind  0.75 


1.1 


Data  for  Computer  Input 

initial  ( TEN  STORY  BUILDING  WITHOUT  VERTICAL  LEAKAGE 


data 

( -20 

20 

bui Iding 

i 10 

1 

heights 

\ 5 

10 

temperature 

f 2 

profi les 

1 

70 

( 2 

1 

50 

outside 

( 2 

1 

pressure 

< 30 

30 

.15 

profi les 

(0.7 

-0.7 

1 


At  standard  conditions  of  21°C 


( 70  °F) 


and  one  atmosphere. 
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1 

2 0 1 

0.75 
0.75 

10  2 

.42 

2.  Example  2 

This  is  a 10  story  building  which  is  70°F  inside.  Outside  the  air  temperature 
is  -5°F  and  there  is  no  wind.  This  building  has  a stairwell  and  an  elevator.  The 
flow  areas  which  are  generally  the  same  vertically  are  listed  in  table  B2  and  the 
flow  exponents  are  taken  to  be  0.5.  The  stairwell  is  pressurized  by  a net  550  cfm 
per  floor.  The  elevator  shaft  has  a 4 ft2  vent  to  the  outside  at  the  top.  On 
floors  2 through  10  the  elevator  lobby  separated  from  the  building  by  doors  that 
automatically  close  in  the  event  of  a fire.  The  flow  coefficient  is  taken  as  0.65 
in  all  connections. 


bui Iding 
data 


10 

1 

0 

1 

2 


10 

0 

.65 

.65 


shaft 

data 


STAIRWELL 
80000  1 

0 550 

1 .65 
0 


Table  B2.  Flow  areas  for  example  2 
Connection  location  Area  (ft2) 


Between 

stairwell 

& building 

0.42 

Between 

building 

& outside 

1.5 

Between 

elevator 

& elevator  lobby 

0.65 

Between 

elevator 

lobby  & building 

0.55 

2.1  Data  for  Computer 

Input 

initial 

(TEN 

STORY 

BUILDING  WITH  ELEVATOR 

data 

(“5 

2 

0 

bui Iding 

fiO 

1 

heights 

1 5 

10 

temperature 

i 1 

profile 

1 1 

1 

70 

outside 

( 1 

1 

pressure 

< o 

1 

1 

data 

1 1 
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bui Iding 
data 


{10 


shaft 

data 


1st  { 

1 

1 

floor  < 

0 

0 

1 

( 

1 

0.65 

0.75 

2 

10 

2 

2nd 

1 

0 

1 

through 

2 

0.65 

0.55 

10th  floors 

1 

0.65 

0.75 

0 

0 

0 

{ 2 

0 


0 

0 


1 


1 


1 


shaft  l 


shaft  2 


STAIRWELL 
80000  1 10 

0 550 

1 .65  .42 

0 


ELEVATOR 


2.7E6 

1 

10 

0 

0 

2 

.65 

.65 

2 

2 

10 

1 

.65 

4.0 

3 

1 

1 

.65 

.65 

1 


1 


2 . 2 Example  2 Output 


The  output 
for  Cases  2 and 


for  example 
3)  is  given 


2 case  1 (the  data  above  not  including  modifications 
in  appendix  C. 
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APPENDIX  C.  EXAMPLE  OUTPUT 
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OUTSIDE  TEMPERATURE 


00 

CN 

O 

o 


uJ 


< 


13 

Q. 


ZD 

O 


UJ 

Q- 


<r 

>c 

UJ 


O 

UJ 

o 


l/D 

UJ 


O 

a 

a. 


UJ 

a 

3 

6- 

< 

a 

UJ 

CL 

s 


oooooooooo 


oooooooooo 


h-  oooooooooo 

X oooooooooo 

OH-  

•-«  u.  in  in  in  m in  to  id  in  m in 
uj  *-oj(r)*^incot^a)(j) 

x 


a»^^iooi»-^ioo>riQ: 

O^NOCONOPHOOUJ 

QMD’J-WOOJMn^H 


oooooooooo 

oooooooooo 
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THE  FOLLOWING  UNITS  ARE  USED  FOR  OUTPUT 


TEN  STORY  BUILDING  WITH  ELEVATOR 

ADJUSTED 

TEMP  FIXED  DIFFERENTIAL  FLOW  FLOW 

FLOOR  COMPARTMENT  PRESSURE  PROFILE  FLOW  CONNECTION  TO  PRESSURE  COEFFICIENT  AREA  FLOW 


h~ 

UJ 

Z 


UJ 

z 


in  n oj  t-  co  co  t t n in  r-  o cn  cm  co  ooo  cn  co  cm  co  cn  o — cj)  in  co  co  cn  — cn  — t***  o m n cv 


r-  t-  « 0(0(0*  oo 

co  cm  in  oo  oo  o coco 

T (O  — CO  T — COCO 


00  cn  CM 
0-0(0 
co  in  — 


co  oo  t o co  • T in 

t t — co  — — — 

co  co  co  in  cm  co  co 


o cn  co  - o — 

id  t id  co  oo 

CM  in  CM  CM  CM 


in'oo 

T (O  CM 

cm  in  co 


cr>  o o — o o 

co  m in  m cm  m 

o co  r-~  o o-  c* 


- o - o o 

in  m m cm  in 

oco  o t- 


- o - o o 

in  in  uo  cm  in 

0(0  o t Is- 


— o - o o 

in  m in  cm  in 

0(0  o T t" 


- o - o c 

in  in  in  cm  in 

oco  o t c"- 


•n  <3-  o 
cn  cn  — 
O -D  - 

— — CM 


in  cl  t 

oo  cn  in  'no 

(TOD  T (0 


in  in 
co  c*  in 

TOO 


T T 
CO  O 
T -3 


t in  lo  t <? 

oo  cn  in  o o cn 

TOO  ^ vO 


t m in  t «r 

oo  cn  m oo  cn 

n O D ST  O 


t-  in  in 
co  cn  in 

T O CD 


CM  (0  (O 
00  CO  O 
- - O 


O N CO  o O 

t"  cn  o in 

o-o  o o 


0)Nh  <n  CM 

in  — o in  t 

o cm  o o o 


oo  in  cm  oo  in 

T co  — T oo 

o cm  o o o 


CD  CM  cn  CO  00 

CO  If;  — CO  CM 

o CM  o o o 


<n  r- 

CM  (D  CM 
O CM  O 


CM  -I 
UJ  UJ 
3 O 
CL  CL  *-* 
O *-*  (0 
O < »- 
-M-  D 
U-  L/D  O 


cr  < 
o > 

o UJ 


o *-« 
-J  o 
CO  -j 

UJ  UJ 
3 O 
cc  ct  •— < 
O CO 
O < H 
-I  H D 

u.  m o 


co  a: 
o 
»- 
cc  < 
o > 

O UJ 


o a 
u -* 
-i  a 
T -* 

UJ  UJ 
3 O 
cc  o:  m 

o *-*  in 
o < »- 

-H-D 

u-  in  o 


□C  < 

o > 

O UJ 


3 O 

o: 

*-*  in 

< t- 
h-  ZD 

m o 


cn  < 

o > 

o ui 


s i 
O i 
O < 
-J  l 
CD  — 1 
UJ  i 

3 l 
o:  ck  * 
O -*  c 
O < » 


CO 

o 

cn 


T 

(O 


CM 

CM 

in 


oo 


<n 


oo 

cn 

o 


CM 


CM 


CM 


CM 


CM 


CM  CO 


CO  T 


t in 


in  co 


(0 
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.069  1 0.  FLOOR  6 COMPARTMENT  1 .029  143-1.  .051  ?4b 

ELEVATOR  -.021  169-1.  .650  -24:3 


TEN  STORY  BUILDING  WITH  ELEVATOR 

ADJUSTED 

TEMP  FIXED  DIFFERENTIAL  FLOW  FLOW 

FLOOR  COMPARTMENT  PRESSURE  PROFILE  FLOW  CONNECTION  TO  PRESSURE  COEFFICIENT  AREA  FLOW 


ID 

2 


ID 

2 


UJ 

2 


UJ 

2 


cm  cm  cm  <t  cm  co  'T  to  cm  co  «t  cm  to  o to  cm  to  co  cm  co  o o co  in  cm 


CO  O CM  • CO  00 

O Q)  t"  O O 

CM  If)  CO  CM  CM 


CO  CO  if)  . CO  00 
to  O)  CM  to  U) 

lO  <T 


'T’TO) 
CM  O P» 
tO  <T 


<r  « in  cm  • min  • 

CM  CM  (''**-  co  t>* 

• to  in  • 


- o o 

Lf)  CM  a") 
O'TN 


— o - o o o 

Lf)  lo  if)  cm  in  m in 

o to  o <r  o to 


- o o 
in  cm  m 
o r- 


- o o o 

in  in  m cm  in 

o to  O^TN 


- o 
in  in 
o to 


*r  if)  if)  rr  *r 

CO  CD  If)  co  O) 

'TOO  to 


«r  in  in 
•*o  o f) 
%*  o J) 


<r  ’T  ld  if) 

co  o co  o',  in 

N to  ^ O CD 


*r  ^ in  in  o- 

CO  O CO  CD  if)  ro  CD 

or  to  or  o cd  to 


«-  - to  «-  in 

CM  CO  CO  CM  *- 

O CM  o o o 


<T  'TO 

»-  CD  <T  *-  -r- 

O CM  o o o 


CO  T O 

o o to 
O CO  o 


co  in 
o o 
o o 


co  co  in 

O ■*“  f** 
O CO  o 


CO  CM 

o o 
o o 


CM 


CM 


CM 


CM 


2 

ID 

S 2 

H-  O 

CL  •— 

< h- 

CL  U 

S ID 

O cl 

(J  — 

-J  Q 
-i 

ID  ID 

3 Q 

CL  CL 

o *-•  t n 

O < •- 

JHD 
u_  lO  O 


z 

ID 


2 

ID 


z 

ID 


2 

2 

z 

2 

H- 

f- 

o 

)- 

CL 

cc 

*— « 

CL 

< 

< 

h- 

< 

CL 

CL 

(J 

CL 

2 

2 

UJ 

2 

O 

O 

CL 

O 

(J 

CJ 

U 

r-  cl 

-J  o 
CD  -J 

co  tr 

o 

h- 
CL  < 

a 

ID  ID 
3 Q 
cr  — « 

O 

h- 

cr  < 

O > 

o 

►-«  l/> 

o > 

o UJ 

o 

< »- 

o UJ 

-J  _J 

_J 

— J _J 

li_  UJ 

lx. 

L/)  o 

U.  UJ 

z 

ID 

5 z 
►-  o 
o:  *-• 

< H- 
Q-  U 
5 ID 

O a 

u 

-i  co 
O)  -J 
UJ  UJ 
3£  Q 

a:  a:  *-• 
o ^ to 
o < h- 

DH  3 

u.  t/>  O 


Z 

UJ 

2 


CL 

< 

CL 

2 

O 

u 


O)  a 

o 


a:  < 
O > 
o UJ 


2 2 

»-  O 

cr 


O a 

o -• 

-l  Q 
O -J 

c-  UJ  ID 

3 a 
a a « 
o if) 

o < »- 

-Jh  D 

u.  D)  O 


O 

CJ 

o a 

tr-  o 
>- 
o:  < 
O > 
O uj 


Is*  (D 


CD  CD 


CD 
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TEMPERATURE  PROFILE  1 

SHAFT  FLOW  COEFFICIENT  80000. 
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c cn  o ^ o O n cn  • 
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Q UJ 
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000000000)0 

oooooooooo 


»-  2 o 
o o 

DJU. 
“3  U.  U. 
O UJ 

< o 

u 


^r*r<T<3,«T<r<T^*T*rco 

ooooooooooo* 

iD.OlDlOlDlOiOlOlDlO^ 

— «-  *-  o 


Z D 
uj  cn 
a cn 

UJ  UJ 

u.  cc 
u_  a 
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z z 
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o o 
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1 1 1 III  III  111  III  III 

ssssss 
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z z 

Ul  Ul 

s s 
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CL  CL 


z z : 
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cc  cc 
< < 
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o o 
o o 
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CL  CL 
2 2 
o o 
o o 
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MAIN  PROGRAM 


£NBS*PLI BS.SHOW  A. MAIN 
C 

C COMPUTER  PROGRAM  FOR  AIR  FLOW 

C SPECIFICALLY  FOR  ANALYSIS  OF 

C 
C 

c 
c 
c 
c 


ANALYSIS  IN  BUILDINGS 
SMOKE  CONTROL  SYSTEMS 


FROGRAM  VARIABLES 


IS  ACCEPTABLE 


OF  FLOOR 


AI  LEAKAGE  AREA  OF  INTERNAL  CONNECTION 

AO  LEAKAGE  AERA  OF  CONNECTION  TO  OUTSIDE 

C FLOW  COEFFICIENT  BETWEEN  BUILDING  POINTS 

CO  FLOW  COEFFICIENT  TO  OUTSIDE 

CS  FLOW  COEFFICIENT  OF  SHAFT 

E LIMIT  WITHIN  WHICH  CONVERGENCE 

F NET  FLOW  INTO  POINT  I 

FC  FLOW  BETWEEN  INTERNAL  FOINTS 

FF  FIXED  FLOW  INTO  POINT  I 

FO  FLOW  TO  OUTSIDE 

FSS  NET  FLOW  INTO  SHAFT  IS 

F HEIGHT  FROM  GROUND  TO  MIDPOINT 

IB  U G OUTPUT  VARIABLE 

ICON  V INTEGER  USED  IN  SUBROUTINES  BLDGP  AND  SHAFTP 

IF  ICONV  = 0 THEN  THE  PRESSURES  WERE  UNCHANGED 
IFLCCR  FLOOR  LEVEL  WHERE  POINT  IS  LOCATED 
IT  POINTER  TO  TEMP  PROFILE  FOR  POINT  I 

ITS  POINTER  TO  TEMPERATURE  PROFILE  OF  SHAFT 

JC  POINT  NO.  CONNECTED  TO  POINT  I 

JOC  DIRECTION  OF  OUTSIDE  CONNECTION 

N NO.  OF  BUILDING  COMPARTMENTS 

NC  NO.  OF  INTERNAL  POINTS  CONNECTED  TO  POINT  I 

NCO  NC.  OF  OUTSIDE  CONNECTIONS 

NFS  1 BOTTOM  FLOOR  OF  SHAFT 

NFS 2 TOP  FLOOR  OF  SHAFT 

NH  NO.  OF  FLOORS 

NPC  NO.  OF  OUTSIDE  PRESSURE  PROFILES 

NS  NO.  OF  SHAFTS 

NS  1 I VALUE  FOR  START  OF  SHAFT 

NS2  I VALUE  FOR  END  OF  SHAFT 

NT  TOTAL  NO.  OF  POINTS  (BLDG  ANO  SHAFT) 

NT P NC.  OF  TEMPERATURE  PROFILES 

P PRESSURE  AT  POINT  I 

PFO  OUTSIDE  PRESSURE  FROFILES 

FO  OUTSIDE  PRESSURE 

PS  PRESSURE  PROFILE  OF  SHAFT  - WORKSPACE 

PZ  PRESSURE  DUE  TO  ELEVATION  DIFFERENCE 

T TEMPERATURE  FRGFILE  ARRAY 

TITLE  PROJECT  TITLE 

TITSH  SHAFT  TITLE 


FROGRAM 

PARAMETERS 

MB 

MAX 

NO. 

OF 

BUILDING  COMPARTMENTS 

MM 

MAX 

NO. 

OF 

POINTS 

MS 

MAX 

NO. 

CF 

SHAFTS 

MC 

MAX 

NO. 

OF 

CONNECTIONS  FOR  ANY  POINT 

MPO 

MAX 

NO. 

OF 

OUTSIDE  PRESSURE  PROFILES 

MTP 

MAX 

NO. 

CF 

TEMPERATURE  PROFILES 

MFL 

MAX 

NO  . 

OF 

FLOORS 
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MAIN  PROGRAM 


PARAMETER  (MM=140  > MS  = 8 , MC  = 9 , MFC=2 . MT P=2 » MFL= 25 * MB=50 ) 

CCMMCN  NT  * P( MM  ) ,C(MM, MC) »NC (MM) , JC ( MM, MC ) , ITS(MS) , 

1 FC(MM.MC) » P Z ( M M , MC ) , PO ( MM , MPO ) ,CO( MM ,MPO) ,F( MM) , PFO ( MF L . MPO ) , 

2 FF ( MM  ) »FO( MM  » MPO)  ,CS( MS) » PS ( MFL ) * NS1 (MS)  ,NS2 (MS) . 

3 FSS (MS)  «N  » NS  * NPO • ICONV, E, I BUG, AI ( MM, MC ) , AO( MM, MPO) , T I T SH (M S , 5 ) , 

4 NF,  KMFL),  IFLOOR  (MM)  ,T(  MTP  , MFL  ) ,NFS1  ( MS)  ,NFS2(  MS)  , IT  ( MB  ) .NTP 

5 ,NCC( MM) , JOC( MM, MPO) .TOUT 
DOUBLE  PRECISION  P.PO.PS 
COMMCN  /RUN/I RUN 

Cl MEN  S ION  Bl( MM,MC) ,B2(MM,MF0) 

N I TER=5000 
IR  UN  = l 

CALL  INPUT  TO  READ  DATA 

CALL  INPUT 

E=  0,2 
1CS=  1 

SAVE  A I ( I , J ) IN  BKI.JI  AND  FIND 
MAX  VALUE  OF  A I ( I , J ) 


AZ  Z = C 
AMAX  = 0 

DO  10  1=1,  NT 
CO  8 J = 1 , MC 
El  ( I « J > = A I ( I , J) 

IF(AI(I,J)  .GT.  A MAX ) AMA  X=AI  ( I • J ) 

CONTINUE 

DO  S J = 1 » MPO 

E2(  I ,J  ) = AO( I, J) 

IF (AC ( I, J ) .GT.  AMAX) AMAX=AC( I , J ) 

CONT  INUE 
0 CONTINUE 

ADJUST  FOR  LARGE  VALUES  OF  FLOW  AREA 

IF (AMAX  .LT.  0.2) GO  TO  25 
AZ  Z = 1 

AM=  0 .2/ ( AMAX— 0 • 1 ) 

EB  = 0 • 1 * ( 1.0-AM) 

CO  15  1=1,  NT 
CO  12  J= 1 » MC 

IF(A1(I«J)  .LT.  0.1)  GO  TO  12 
AI ( I . J ) = AM*AI ( I « J ) +BB 
2 CONTINUE 

DO  14  J= 1 , MPO 

IF(AC(I,J)  .LT.  0.1)  GO  TO  14 
AO  < I , J ) = AM*AO(  I . J ) +BB 

4 CONTINUE 

5 CONTINUE 

TEMPERATURE  CORRECTION 
25  CALL  CCRR 
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MAIN  PROGRAM 


C 

C CALL  INIT  TO  INITIALIZE  PRESSURE  ARRAY  , p 

C 

22  CALL  INIT 

C 

c 

C DO  LOOP  TO  30  IS  ITERATIVE  SOLUTION  TO  PRESSURE  ARRAY 

C 

24  DO  3C  ITER=1. NITER 

C 

C CALL  BLDGP  TO  SOLVE  FOR  BUILDING  PRESSURES 

C 

CALL  BLDGP 
IC  E = ICON V 

IF ( I CB  .EQ.  0 .AND.  ICS  .EQ.  0 ) GO  TO  40 
C 

C CALL  SHAFTP  TO  SOLVE  FOR  SHAFT  PRESSURES 

C 

CALL  SHAFTP 
ICS=  1C0NV 

IF ( I C B .EQ.  0 .AND.  ICS  .EQ.  0 ) GO  TO  40 
C 

C CALL  PZAD  TO  CALCULATE  PZ  TERMS 

C 

CALL  PZAD 
30  CONTINUE 

C 

C IF  ROUTINE  FAILS  TO  CONVERGE  IN  NITER 

C ITERATIONS  PRINT  ERROR  MESSAGE 

C 

HR  I TE ( 6.800) 

AO  CONTINUE 

HRITE(6»801)ITER 

IF ( A 2 Z .EQ.  0 . ) GO  TO  42 

AZZ  = 0. 

DO  6C  1=1. NT 
DO  5 C J=1.MC 
50  AI  ( I , J )=  B 1 ( I,  J ) 

DO  55  J=1 , MPO 

55  AO  ( I , J ) = B2(  I, J ) 

60  CONTINUE 

CALL  CCRR 
GO  TC  24 
C 

c 

c 

C CALL  OUT  TO  OUTPUT  SOLUTION 

C 

A2  CALL  OUT 

C 

wRiTE(e.eos) 

STOF 

c 

c 

C FORMAT  STATEMENTS 

C 

£00  FORMAT (/////5X. 35 < 1H1  )//5X. 

+35HF  A I LURE  OF  MAIN  PROGRAM  TO  CONVERGE  / /5X , 35 ( 1 Hi  )// ) 
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MAIN  PROGRAM 


801  FORMAT!  I OX , 1 5. 5X , 1 1 H I 7ERAT I C NS  ) 

£05  FORMAT  (11-1) 

END 

aFDG.P  SUBROUTINE  INPUT. L.l 
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SUBROUTINE  INPUT 


$NBS*PLI B$ .ShGW  A. INPUT 
SUBROUTINE  INPUT 

THIS  ROUTINE  READS  AND  PRINTS  DATA 
ANC  INITIALIZES  PZ  ARRAY 

PARAMETER  ( MM  = 14  0 . MS  = 8 . MC=9 . MFO=2 , MT P=2 , M FL= 25 , M B=50 ) 

COMMCN  /PZZ/  PGZ 

COMMON  / IO/T  ITLE  ( 18)  , I OUT,  IUNIT,NC0MP<  MFL)  ,SNCOMP(  MFL  ) 

COMMCN  NT,  P ( MM ) ,C(MM , MC ) ,NC (MM) , JC ( MM, MC)  , ITS( MS)  , 

1 F C ( M M *MC)  , PZ (MM, MC),PC(  MM, MFC), CO ( MM , MPO ),F(MM), PFO( MFL, MPO ) , 

2 FF ( MM  ) , FO ( MM , MPO  ) ,CS( MS) , PS( MFL)  ,NS1 ( MS)  ,NS2 ( MS ) , 

3 FSS(MS),N,NS,NPO,ICONV,E.IBUG,AI(MM,MC),AO(MM,MPO),TITSH(MS,5), 

4 NH, h( MFL) , IFLOQR ( MM) , T ( MT P , MFL ) , NFS 1 ( MS) , NFS 2 (MS), IT(MB)»NTP 

5 , NCC (MM  ) , JOC ( MM, MPO) , TOUT 
COLELE  PRECISION  P.PO.PS 
CHARACTER  PAR  *6 

DIMENSION  II( MFL)  ,TT( MFL) , PAR  < 7 ) , CW ( MPO ) , PH( MFL) ,NZZ(MM ) 

CATA  FAR/'  MM','  MS*,'  MC • , • MFO * , • MTP • , • MFL* , • MB'/ 

IBUG=0 

READ  AND  WRITE  PROJECT  TITLE 

READ (5,600) (T ITLE ( I ) . 1 = 1 . 18) 

WRITE(C,e01)( TITLE ( I)  ,1  = 1 ,18) 


READ  GENERAL  DATA 


TOUT  = OUTSIDE  TEMPERATURE 
IUN  IT  = 1 FOR  SI  UNITS 
= 2 FOR  ENG  UNITS 
ICLT  = 0 FOR  NO  SUMMARY  OUTPUT 
OThERWlSE  IOUT  IS  FILE  NO.  TO 
WhICH  SUMMARY  OUTPUT  IS  WRITTEN 

RE AC <5,7 00) TOUT, I UNIT, I OUT 
Wrt  ITE(€,41 1 )TOUT»  IUNIT.IOUT 

IF ( ILNIT  .GT.  2 .OR.  IUNIT  .LT.  1 ) GO  TO  105 
READ  hEIGHTS 

N N = 0 FOR  INPUT  OF  ALL  hEIGHTS 
C NN  = 1 FOR  CALCULATION  OF  hEIGHTS 

C 

READ (5 ,700) NH , NN 
WR  ITE (C,412)NH,NN 
I F ( N h »LE.  MFL  ) GO  TO  89 
IP AR  =6 
GO  TO  110 

69  IF (NN  .EG.  1 ) GO  TO  97 

READ (5 ,700) (H( I ), I =1 ,NH) 

WR  ITE(6«413)(H(I)  , 1 = 1 , NH ) 

GO  TO  99 

97  READ  (5 ,700  )H(  1 ) ,DH 
wR  I T E ( e,414)H(  1 ),DH 
CO  9 6 1=  2 , NH 
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SUBROUTINE  INPUT 


I M = I—  1 

58  MI)  = MIM)+DH 
READ  TEMPERATURE  PROFILES 

59  REAC (5»700)NTP 
WRITEC  6.415JNTP 

IF ( N T P .LE.  MTP ) GO  TO  SO 

IPAR=5 

GO  TC  110 

50  DO  2 IF=1,NTP 

REAC (5  ,700)NNN, (I  I ( J ) , TT  ( J ) , J = 1 ,NNN) 

WR  ITE(6»416)NNN»(  I I(J  ) »TT(J) , J = 1 ,NNN) 

IF ( NHN  .GT.  1 ) GO  TC  2 
CO  1 I FF  = 1 » NH 

1 T(  IP .IFF )=TT(  1 ) 

GO  TC  3 

2 J=  1 
JP  1 = 2 

CO  4 I FF  = 1 9 NH 

T(  IP  , IFF  ) =TT ( J ) + ( TT ( JP 1 ) — TT ( J ) )*( IFF- I I( J) )/( II ( JPI  )-I I ( J)) 

IF  ( I FF  ,NE.  I I ( JP I ) ) GC  TC  4 

IF ( JF  1 .EQ.  NNN )G  0 TO  4 

J=JP  1 

JPl  = J+l 

4 CONT  1NUE 

2 CONTINUE 

C 
C 

C READ  OUTSIDE  PRESSURE  PROFILES 

C NN=0  FOR  INPUT  OF  ALL  PRESSURES 

C NN=1  FOR  CALCULATION  BY  POWER  LAW 

C 

REAC (5.700  JNPO.NN 

WR  XTE( 6»417)NPC  »NN 

IF ( N FO  .LE.  MPO )GO  TO  SI 

IP AR  = 4 

GO  TC  110 

51  IF ( N h .EQ.  1 ) GO  TO  81 

C 

C READ  ALL  OUTSIDE  PRESSURES 

C 

CO  6 1=1. NPO 

€ READ (5,7  00 ) PGZ . ( PFO ( J . I ) , J= 1 , NH ) 

WR ITE(€.413)PGZ,( PFO ( J . I ) » J= 1 » NH ) 

GC  TC  es 
C 

C CALCULATE  OUTSIDE  PRESSURES 

C PATMCS  IS  ATMOSPHERIC  PRESSURE  (PA) 

C 

8 1 REAC(5.7CO)VW,hW,XW,(CW< I ) .1  = 1 ,NPO) 

WRITE( e.419)VW,HW,XW,(CW(I) *1=1, NPO ) 

IF ( IUN  IT  .EQ.  1 )VW  = VW*0 .2778 
IF ( I UN  IT  .EQ.  2)VW  = VW*0. 4470 
PATMCS=10 1325. 

TO  0=1CLT  + 273. 

IF ( IUN  IT  .EQ.  2 )TOO=( TOUT  +460 . ) /I .6 
PVA=  17  8. 4*VW*VW/T0C 
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SUBROUTINE  INPUT 


Z=-0  .034  17/TOO 

IF ( I UN  IT  .EQ.  2)Z=0.3048*Z 

CWM=CW  ( 1 ) 

IF ( N FQ  .EQ.  1 ) GO  TO  212 

DO  211  1 = 1. NP  O 

IF(CMI)  .LT.  CWM)CWM=CW(  I) 

211  CONTINUE 

212  PG  Z=  F ATMOS* EXP  < H(  NH  ) *Z  ) +CW  M*PV  A*  C ( H (NH  ) /HW  )**  ( 2 . *X  W ) )-100. 
CO  210  1=1. NH 

PH ( I )=PATMOS*EXP( H ( I ) 4Z ) 

210  CONTINUE 

DO  82  1=1. NPO 

CO  8 2 J=  l . NH 

PFO ( J * I ) =PH( J) +CW(  I ) *PVA*( ( H ( J )/HW)** (2 .*XW ) )-PGZ 
82  CONTINUE 

C 

c 

C BUILDING  DATA  INPUT 

C NFLS  = NO.  OF  FLOORS  IN  BUILDING 

C IF  1 = LOWER  FLOOR  IN  SERIES  OF  SIMILAR  FLOORS 

C I F 2 = UPPER  FLOOR  IN  SERIES  OF  SIMILAR  FLOORS 

C NCC  = NO.  OF  COMPARTMENTS  PER  FLOOR 

C NZ  = NO.  OF  CONNECTIONS  TO  COMPARTMENTS  ON  SAME  FLOOR 

C NA  = NO.  OF  CONNECTIONS  TO  COMPARTMENTS  ON  FLOOR  ABOVE 

C 

85  1 = 0 

SNCCMP  ( 1 ) =0  * 

READ (5  .700) NF L S 
WRITE(6.420)NFLS 
IF ( NFL  S .GT.  NH )GO  TO  106 
7 REAC  (5.700)  IF1  , IF2.N0C 

WR  I T E ( 6 « 40 0 ) I F 1 .IF2.N0C 
IF ( I f 1 .GT.  IF  2 )GO  TO  107 
NCOM  F ( IF  1 ) =NOC 
IFP=  IF  1+1 

SNCO  M P ( I FP  ) =SNC0MP ( I F 1 >+NOC 
CO  1C  IZ=1 . NOC 
1=  1+  1 

READ (5  .700) NZ . NA. NN0.FF( I) .IT  < I) 

WR IT E( 8.401 )NZ,NA .NNO.FF (I ) .IT (I ) 

NZ Z ( I ) =NZ 

NN=N  2 + NA 

IFLOCR ( I ) =1 F 1 

IF (NN  .LE.  MC)G0  TO  111 

IP  AR  = 3 

GO  TC  110 

111  IFCNNC  .LE.  MPO )G  0 TO  112 
IP AR  =4 

GO  TC  110 

112  IF ( IT ( I ) .GT.  NTP  .OR.  IT(I)  .LT.  1 ) GO  TO  102 
NC ( I ) = NN 

IF  (NZ  .EQ.  0)  GO  TO  83 

INFUT  CONNECTIONS  TO  COMPARTMENTS  ON  SAME  FLOOR 

REAC  (5.7C0) ( JC( I. J ) .C< I * J) .AI< I . J)  . J=1 .NZ) 

WRITE (6 .402) 

WRITE(6.403)(JC(I.J).C(I.J).AI(I,J).J=1.NZ) 
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SUBROUTINE  INPUT 


CO  6 2 J=I.NZ 

62  JC< I . J )= JC( I. J)+SNCCMP< IF1 ) 

63  IF  <NA  .EC.  0)  GO  TO  6 
C 

C INFUT  CONNECTIONS  TO  COMPARTMENTS  ON  FLOOR  ABOVE 

C 

NP  = N 2+1 

READ (5 ,700 H JC ( I , J ) ,C ( I . J ) . AI ( I . J ) , J= NP, NN) 

WR ITE( 6,404) 

WRITE (6,403) ( JC<I,J),C(I.J),AI(I.J) . J=NP, NN) 

DC  66  J=NP,NN 

6 6 JC ( I , J ) = JC( I , J ) +NCCMP ( I F 1 ) + SNCO MP ( IF1  ) 
e NCC ( I ) =NNO 

I F ( N NO  .EQ.  0 ) GO  TO  10 
C 

C INFUT  CONNECTION  TO  OUTSIDE 

C 

REAC  <5,700)  (JOC(I * J J ) »CQ(I»JJ)  «AO( I * J J ) , J J = 1 , NNO) 
WR  I T E ( 6 , 40  5 ) 

WR  ITE< 6,40  3 ) ( JCC(  I , J J ) „ C C( I , J J ) , AO ( I , J J ) , JJ  = 1 ,NNO) 
CO  9 J J= 1 * NNO 
J= JOC(  I . JJ) 

9 FO ( I , J J ) =PFO(  IF  1, J ) 

10  CONTINUE 

I F ( I F 1 .NE.  IF 2 )GO  TO  11 
IF  ( I F 1 .EQ.  NFLS) GO  TO  20 
GO  TC  19 
C 

C ASIGN  CATA  FOR  FLOORS  SIMILAR  TO  FLOOR  IF1 

C 

11  IFP=IF1  + 1 

CO  17  IFF=IFP,IF2 
NCOMF  ( IFF  )=NOC 
IFFP= IFF+1 

SNCOMP! I FFP ) =SNCO  MP ( IFF ) +NOC 
CO  16  IZ=1.NQC 
1=  1+  1 

I1  = I2  + SN  CO  M P ( I F 1 ) 

IFLCCR  ( I )=  IFF 
FF ( I )=FF( II) 

IT< I )=  I T < I 1 ) 

NN=NC<  1 1 ) 

NN0=NCC<  ID 
NC ( I ) = NN 
NCC ( I ) =NNO 

IF ( IFF  .NE.  NFLS ) GO  TO  23 
NN  = N 2 Z ( I 1) 

NC ( I ) = NN 

23  IFINN  .EG.  0 ) GC  TO  14 

DO  12  J= 1 , NN 
C< I.D=C<I1,J) 

A I (I , J ) = A I ( 1 1 , J ) 

JC ( I , J )= JC( II  , J )+  SNCOMP { IFF ) — SNCOMP ( I F 1 ) 

12  CONTINUE 

14  IF (NNO  .EQ.  0 ) GO  TO  16 

DO  16  JJ=l,NNO 
JOC ( 1 , J J ) = JOC< 1 1. J J) 

J= JCC(  I , J J) 


34 


SUBROUTINE  INPUT 


CO ( I ,JJ  )=CO( II  , JJ  ) 

AO  ( I « J J ) =AQ( I 1 . JJ  ) 

15  FO ( I ' J J )=PFO(  IFF, J ) 

16  CONTINUE 

17  CONTINUE 

18  IF  ( I F 2 .EQ.  NFLSIGC  TO  20 

19  CONTINUE 
GO  TC  7 

20  N=  I 
N2  = N 

IF (N  .LE  . MB) GO  TO  114 
IP AR=7 
GO  TC  110 
C 

C SF AFT  CATA  INPUT 

C 

114  READ (5 .700 )NS 

IF ( N S .LE.  MS ) GC  TC  113 
IP AR  =2 
GO  TC  110 

113  DO  ICO  I S= 1 » NS 

READ  <5.603 ) <T I T SH < I S . I ) • 1= 1 . 5 ) 

WRITE(6,406)( TITSH{ IS. I ) .1=1 .5) 

READ  (5.7  0 0>CS< IS) .NFS1 < I S ) .NF S2( I S ) , I T S< I S ) 

WR I TE ( 6 • 407 ) CS(  IS) » NFS1 ( IS ) ,NF S2 ( IS ) . ITS(  IS) 

N 1 =N  2 + 1 

N2  = N 1 + NF S2 ( IS) -NFS  1C  IS) 

NS  1 ( IS  ) = N1 
NS2 ( IS ) =N2 
IFF  = NFS1 C IS)-1 

READ (5,700) NNO . FF  F.JCP.CC.AA 
WR  ITE(6.408)NNC.FFF, JCP.CC  »AA 
IF ( NNO  .EQ.  0 ) GC  TO  21 

RE AC  (5.7  CO) (JOCCN 1 .J) ,C0(N1 ,J)  .AO(Nl.J)  ,J  = 1 ,NNO) 
WRITE (6 ,403) ( J0C( N 1 , J ) . CO ( N1 . J ) , AO ( N1 , J ) , J= 1 , NNO) 

21  DO  24  I=N1.N2 

NC(  I )=1 

NCO ( I ) =NNO 

IF  F = 1FF+  1 

IFLOCRC I )=IFF 

I F ( IFF  .GT.  NFLS) GO  TO  25 

FF ( I )=FFF 

IF ( J CP  .GT.  NCOMP ( IFF) ) GO  TO  25 
JC ( I . 1 ) = JCP+SNCOM P ( IFF) 

C( I, 1 )=CC 
AI  (I .1  ) =AA 

26  IF ( NNO  .EQ.  0 ) GO  TC  24 

CO  22  J= 1 , NNO 
JJ= JCC(N1 , J) 

FO ( I , J ) = PFO( IFF , J J ) 

JOC<  I , J)  = J J 

CO  ( I .J  ) = C0(N1 , J) 

22  AO (I  ,J  )=AO(Nl , J) 

GO  TC  24 

25  NC(I)=0 

GO  TC  26 
24  CONTINUE 

C 
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SUBROUTINE  INPUT 


EXCEPTIONS  TO  GENERAL  SHAFT  INPUT 
NNN  = NO.  OF  EXCEPTIONS 
KE  = 1 FOR  FF  EXCEPTION 
C KE  = 2 FOR  OUTSIDE  CONNECTION 

C KE  = 3 FOR  INTERNAL  CONNECTION 

C 

READ {5  .7  00) NNN 

IF (NNN  .EQ.  0 ) GO  TO  100 

DO  69  IK=1,NNN 

READ  <5  .700 )KE.  IFF 

MR  I T E ( 6 • AO 9 ) KE  t IFF 

I=NS  1 ( IS )+ IFF-NFS 1 ( IS) 

IF  (KE  .EQ.  1)  GO  TO  41 
IF  ( K E .EQ.  2 ) GO  TO  42 
I F ( K E .EG.  3) GO  TO  51 
C-0  TC  10  4 

41  READ (5  .700)FF ( I ) 

WR ITE ( 6. 410 )FF ( I) 

GO  TC  69 

42  READ (5»700)J . CCG* A AO 
MR  IT  E ( 6.  405) 

WRITE1 6 ,403) J .CCO.AAO 
NN  C=  N CO  ( I) 

IF ( NNC  .EQ.  0 ) GO  TQ  44 
CO  43  K=1,NNC 

IF(JCC(I.K)  .EQ.  J ) GO  TO  46 

43  CONTINUE 

44  NJC=NNC+ 1 
NCO ( I ) =N  JG 

47  PO( I ,N JO )=PFO( IFF , J) 

JOC ( I .NJO)=J 
CO( I .NJO)=CCO 
AO ( I «N  JO )=AAO 
GO  TC  69 
46  NJO  -K 

KK  = K4  1 

IF (CCC  .NE.  0 ) GO  TC  47 
NJO  =N  NC—  1 
NCC ( 1 )=N JO 

IF ( NJO  .EQ.  0 ) GO  TO  69 
DO  49  K=  KK . NNC 
KM=K— 1 

FO ( I ,KM)=PO ( I . K ) 

JOC ( I « KM ) = JOC ( I , K ) 

CO (I  .KM)=CO(I.K) 

49  AO  ( I iKM  =AO(  I . K ) 

GO  TC  69 

51  READ (5 .700 ) JCP .CC . AA 
WRITE  (6. 402) 

WRITE( 6 ,403) JCP.CC .AA 
J= JCF+SNCOMP< IFF) 

NN  = N C ( I ) 

IFCNN  .EG.  0 ) GO  TO  53 
DO  52  K=1,NN 

IF(JC(  l.K)  .EQ.  J ) GO  TO  55 

52  CONTINUE 

IF (CC  .NE.  0. >GC  TC  53 
MR  ITE ( 6. 520  ) IS  .KE . IFF 
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S3 

GO  TC  69 
NJ=NN+1 

£4 

NC ( I 1 =N J 
JC  ( I ,N J )=J 

55 

C< I .NJ )=CC 
A I ( I »N  J ) =A A 
GO  TC  69 
N J =K 

KK=K4 1 

IF  ( A A .NE.  O.iGO  TO 

54 

N J = N N—  1 
NC ( I )=NJ 

IF  ( N J .EQ.  0)  GO  TO 

69 

61 

CO  61  K=KK,NN 
KM  =K—  1 

JC(  I . K M ) = JC  ( I . K ) 
C ( I .KM ) = C< I . K 1 
AI  ( I ,KM ) =A I ( I . K ) 

69 

CONTINUE 

100 

CONT  1NUE 

NT  =N  2 

IF  (NT  .LE.  MM  ) GO  TO 

160 

c 

c 

IP AR  = 1 
GO  TC  110 

PRINT  OUTSIDE  TEMPERATURE 

C 

160  WR  ITE<6.  601  M TITLE!  I)  .1  = 1 . 12) 

IF ( I ON  IT  .EQ.  I )WRITE(6,800)TOUT 
IF ( ION  IT  .EQ.  21WRITE! 6.5C01TOUT 
IF ( I UN  IT  .EQ.  2 )T0UT  = { TOUT-32 . )/l .8 
TOUT  =TCUT  +273. 

C 

C PRINT  FEIGHT  AND  TEMFERATURE  PROFILES 

C 

IF  < I UN  IT  .EQ.  1 )WR ITEC 6, 81 1 ) ( IP. IP=1 .NTP) 

IF ( I UN  IT  .EQ.  2)WRITE(6.511  )(IP.IP=1»NT°) 

HR  I T E ( 6 , 8 1 3 1 
CO  3 C 1FF=  1 »NH 

20  HR ITE( 6.8121H! I FF  ) «(T(  IF»I FF ) * I P= 1 « NT P 1 

C 

C CONVERT  TEMPERATURES  TO  DEG  K 

C 

DO  32  IFF=1.NH 
CO  32  IF  = 1 » NT  P 

IF ( I UN  IT  .EQ.  2)T(  IP, IFF)=(T(  I P . IFF )- 32  . ) /I  .8 
33  T< IP , IFF  ) = T(  IP.  IFF  > + 273 . 

C 

C PRINT  OUTSIDE  PRESSURE  PROFILES 

C 

IF  I I UN  IT  .EQ.  1 ) GO  TO  79 
HRITE<6,514)(  IP  » I P=  1 * NPC  ) 

HR  I T E ( 6 • 8 1 3 ) 

CO  7 € IFF=  1 »NH 
DO  77  J= 1 » NPO 

77  PFO ( IFF. J )=PFO(  IFF , J )/248. 8 

HR ITE(6.515)H( IFF  ) »(PF0C IFF  » J ) »J  = 1 . NP  0 ) 

CO  76  J= 1 . NPO 
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7e  PFC<  IFF, J)=PFO(  IFF,J)*24E.8 

76  CONT  IN UE 

C-0  7C  6 3 

79  WR  ITE(6.814)(  IP.IP=l,NPO) 

WR  I T E ( 6 , 6 1 3 ) 

DO  31  IFF=I ,NH 

WR  ITE(C»  615)H( IFF ) , (PFO< IFF ,J ) . J=1 ,NPO) 

31  CONTINUE 

C 

C CORRECT  FOR  CONNECTIONS  ONLY  INPUTED  ONCE 

C 

63  DO  6C  1=1, NT 

NN=NC<  I) 

IF  { N K .EC.  0 ) GO  TO  60 
CO  56  J J = 1 , NN 
J = JC ( I , J J) 

IF  ( J .EQ.  0 ) GO  TO  58 
NN  J = NC ( J ) 

IF < N N J .EQ.  0 ) GO  TO  57 

DO  56  IA  =1 , NN J 

IF(JC(J,  IA  ) .EQ.  1)60  TO  58 

56  CONTINUE 

57  NNJsl'NJ+l 

IF ( NN  J .LE.  MC)GO  TO  5S 
IPAR  = 3 
GO  TO  110 
£9  NC(J)  = NNJ 

JCC  J ,NNJ )=I 
C ( J , NN  J ) =C  ( I , J J ) 

AI  ( J ,NNJ  ) = AI ( I , JJ ) 

IF  ( J .GT.  N .OR.  I .GT.  N ) GO  TO  58 
PZ ( J *NN J ) =— PZ ( I , J J ) 

58  CONTINUE 

60  CONTINUE 

C 

C CORRECT  UNITS 

C 

IF ( I UN  I T .EQ.  2)CALL  UNITS 

C 

C INITIALIZE  PZ  FOR  BUILD  COMPARTMENTS 

C 

67  CO  4 C 1=1,  N 

NN  = NC ( I) 

IF  (NN  .EG.  0 } GO  TO  40 
IA  = I 1(1) 

IF  1=  IFLOOR  ( I) 

CO  36  J J = 1 , NN 
J=JC  ( I ,J J) 

IF J= I FLOOR! J) 

IF ( 1 F I .EQ.  IF J )GQ  TO  36 
IB=IT( J) 

TEMPA=0.5*(T(  IA,IFI)+T<  IB.IFJ)  ) 

PZ ( I , J J) =3462. *(H(  IF J)-H(IFI  ) ) /TEMP A 
38  CONTINUE 

40  CONTINUE 

C 

C INITIALIZE  PZ  FOR  SHAFTS 

C 
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DC  SC  IS  = 1 »NS 
N1=N£1 (IS) 

N2=NS2( IS)-I 
IT T=  IT S ( IS) 

DO  4 S I=N1,N2 
IF  1=  IFLOOR  ( I ) 

IF J= IF  1+  1 

TEMPA=0.5* (T(  ITT,  IFI)+T( ITT,  IF  J)  ) 

PZ ( I . 1 )=2462.*(  H<  I FJ)-H( IFI ) ) /TEMP A 
45  CONTINUE 

50  CONTINUE 

c 

C CHECK  SHAFT  CONNECTIONS 

C 

CO  240  IS=1.NS 
N1  = N S 1 (IS) 

N2=NS2( I S) 

CO  229  I =N  1 . N2 
NN=NC(  1 ) 

IF  ( N F .EQ.  0 ) GO  TO  239 
CO  226  J=1.NN 
JJ  = JC(  I, J) 

IF ( I FLOOR ( I ) .NE.  IFLOOR ( JJ )) GO  TO  103 
236  CONTINUE 

239  CONTINUE 

240  CONTINUE 
RETURN 

C 

C 

C DIAGNOSTIC  OUTPUT 

C 

102  WRI 1E(  6.902 ) I . IT(  I ) 

GO  TO  109 

103  KR ITE( 6 .903) 

GO  TC  109 

104  taRITE(6.904) 

GO  TC  109 

105  \»R  ITE(6*905) 

GO  TC  109 

106  1*R  ITE(6.906) 

GO  TC  109 

1C7  1*R  ITE(  6.907) 

GO  TC  109 

110  \#RI  TE(  6.910  )PAR(  I PAR) 

C 

C PRINT  CORRECTED  BUILDING  DATA 

C 

109  *RI 7E( 6.940) 

DO  7 C 1=1.  N 
NN  =NC ( I) 

IF ( N N .GT.  0 ) GO  TO  180 

*RITE( 6. 941 ) I . IFLOOR ( I ) , IT ( I ) , FF ( I ) 

GO  TC  182 

180  *R ITE( 6.942) I , I FLOOR! I ) . I T( I ) , F F ( I ) , JC( I . 1 > ,C ( I . 1 > . A I < I , 1 ) 
IF ( N F .EG.  1 ) GO  TO  182 

V*R  ITE  (6  »942)(  JC(  I , J ),C  (I  . J ) .AI  ( I . J)  .J  = 2 ,NN) 

182  NNC  =FCC ( I ) 

IF ( NFO  .EQ.  0 ) GO  TC  70 
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70 

C 

C 

c 

72 

74 

75 

£0 

C 

C 

C 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 

413 

4 14 

415 

4 16 

417 

418 

419 

420 

500 

£11 

£14 

£15 
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WRITE(6,944){ JOC(  I , J ) , CO ( I , J ) , AO ( I , J ) , J=1  ,NN0) 
CCNT  INUE 


PRINT  CORRECTED  SI-AFT  INPUT  DATA 


CO 

8 C 

IS  = 1 , 

NS 

kR 

ITE 

( 

6,816 

) ( TITSH 

(IS, 

I ) 

• 

I 

L ,5 

) 

Vi/R 

ITE 

( 

6.806 

)IS.CS( 

IS), 

IT 

c 

( 

I : 

S) 

N1 

=N£ 

1 

( I S) 

N2 

= NS 

2 

(IS) 

HR 

ITE 

( 

6,807 

) 

CO 

7 £ 

M 

II 

z 

« 

N2 

NN 

= NC 

( 

I) 

IF 

(NP 

.GT. 

0 ) GO  TO 

72 

HR 

ITE 

( 

6 ,801 

) IFLCGR 

( I)  , 

FF 

( 

I 

) 

GO 

TC 

74 

HR 

ITE 

( 

6,808 

) IFLCOR 

( I ) , 

FF 

( 

I 

) ■ 

. JC 

(I 

, 1) 

,C(  I 

.1  ) , A I (1,1  ) 

IF 

(NN 

.EG. 

1 ) GO  TO 

74 

HR 

ITE 

( 

6 . 809 

)(  JC(I  , 

J)  ,C 

( I 

* 

J 

) 

i A I 

(I 

, J) 

, J=2 

.NN) 

NNC  = NCC(  I ) 

IF  (NNO  .EQ.  0 ) GO  TO  75 

HRITE(6,810)( JOC( I . J ) . CC ( I , J ) . AO ( I . J ) , J=1  * NNO ) 

CONT  1NLE 
CONT  INUE 
STOP 

FORMAT  STATEMENTS 

FORM  AT  (5X,  5HIF1  =.  13,71-1,  IF2  =,I3,7H,  NOC  =,I3) 

FORMAT (5X  » 4HNZ  =,I3,6H  NA  = ,I3,7H,  NNO  = .I3.6H,  FF  =,F8«1, 

+ 7 P , IT  = , 13) 

FORM  AT (5X, 25HCONNECTION  ON  SAME  FLOOR  ) 

FORMAT  (5X.3HJ  =,I3.5H,  C =.F10.3,5H,  A =,F9.4) 

FORM  AT (5X, 26HCGNNECT ION  TO  FLOOR  ABOVE  ) 

FORM  AT  <5X, 22HC0NN ECT ION  TO  OUTSIDE  ) 

FORMAT (5X.5A4) 

FORMAT (5X.4HCS  =.F9.l,8H,  NFS1  =,I3,8H,  NFS2  =,I3,7H,  ITS=,I3) 
FORMAT  <5X, 5HNN0  =,I3,7H,  FFF  =,F8.1,5H,  J =,I3,5H,  C =,  F10.3, 

+ £ H , A = , F9 • 4 ) 

FORMAT  (5X.4PK.E  =,I3.  7Pi,  IFF  = ,I3) 

FORMAT (5X  ,4HFF  - , F 8 • 1 ) 

FORMAT (5X.6PT0UT  =,F6.C,9H,  IUNIT  =,I3.8H,  IOUT  =,I3) 

FORMAT (5X,4  HN  H =,I3,6H,  NN  =,I3) 

FORMAT (5X »7HHE IGHT S /(10F8.2)) 

FORMAT  (5X,6HH(  1 ) =,F8.2,6H,  DH  =,F8.2) 

FORM  AT  (6X, 5HNTP  =,I3) 

FORMAT (5X.20HTEMPERATURE  PROFILE  / 15 , ( 1 0 ( 14 , F7 . 1 ) ) ) 

FORMAT  (5X.5PNPO  =,I3,6h,  NN  =,I3) 

FORMAT (5X,5HPGZ  = , F12 . 1 /17HPRESSURE  PROFILE  /(10F12.1)) 

FORMAT (5X  ,4HVW  =,F6.1,6H,  HW  = ,F6.1,6H,  XH  =,F4.2,6H,  CW  =, 

+ ( 10  F4  .2  ) ) 

FORMAT </5X, 6HNFLS  =,I3) 

FORMAT (//10X,20HOUTSIDE  TEMPERATURE  .F6.1.2H  FI 

FORMAT!  ///5X  ,6liFEICHT  ,5X»  29  H TEMPERATURE  PROFILES  (DEG  F)  / 

+ 7X, 2FFT ,3X, 1 916) 

FORMAT (////5X .6HHEIGHT  , 5X , 26H0UTSI DE  PRESSURE  PROFILES 
1 IIP  (IN  H20  ) / 7X , 2P.F  T«3X«8I10) 

FORMAT (Fll ,2,3X,8F10«3 ) 
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520  FORMAT  !///5X»  15FERRGR  IN  SHAFT  , I 2 . I 5HEXCEPT  I ON  KE  - ,12, 

♦ 2 X » EH FLOOR  ,13//) 

£00  FORMAT ( 1 8A4 ) 

£01  FOR M AT ( l F 1/// 1 OX » 18A4///) 

£03  FORMAT (5A4) 

700  FORMAT ( ) 

£00  FORM  AT (//I  OX  » 20HOU  T SIDE  TEMPERATURE  ,F6.1,2H  C) 

£01  F0RMAT(I13iFll.l) 

£06  FORMAT  < 10X  * 1 2HSFAFT  NUMBER  , I 4/ 1 OX , 1 7H SHAFT  COEFFICIENT  .F10.2/ 

1 10X  .20HTEMPERATURE  PROFILE  .14) 

007  F0RMATC/21X.5HFIXED.25X. 4HFL0W  * 12X . 4HFL0W / 1 OX , 5HFL00R , 6 X , 

1 4FFL0Vy,5X,12HCGNNECTED  TO  *6X  . 1 1 HCOEFFIC  IE  NT  .6X.6H  AREA 

2 / ) 

e08  FORMAT ( I 13, FI  1.1, 6X ,5HPG INT , I 5 ,F 16 . 1 , F 1 5. 4) 

ECS  FORMAT ( 3 OX , 5HP0 IN T . I 5 . F 1 6 . 1 .FI  5. 4) 

eiO  FORMAT (30X .7H0UTS IDE  , I 3 , FI 6 . 1 . F 1 5 . 4 ) 

£11  FORM  AT ( ///5X.6FFEIGHT . £ X , 29HTEMPERA TURE  PROFILES  ( DEG  C)  / 

+ 7 X , 2HM  ,3X, 1 9 16) 

£12  F0RMATCFH.2.3X.19F6.  1 ) 

£13  FORMAT!/) 

£14  FORMAT (////5X .6HHEIGHT  . 5X . 26H0UTSI DE  PRESSURE  PROFILES 
1 1 1 F (PASCALS)  /7X.2HM  .3X.8I10) 

£15  FORM  AT ( F 1 1 • 2 . 3X ,8  F 1 0 . 1 ) 

£16  FORMAT <///10X . 5A4 ) 

817  FORMAT! IOX.45HFLOW  COEFFICIENTS  CORRECTED  FOR  TEMPREATURE  ) 

502  FORMAT (10(/)«10X,  1 1HCC  MPARTMENT  ,14/ 

1 1 OX  ,2 OHTEMPERATURE  PROFILE  .I4.17H  DOES  NOT  EXIST  / 

+ 1 OX  , 1 6 H PRO GRAM  STOPPED  .10!/)) 

503  FORMAT (10!/) ,5X .23HSHAFT  CONNECTION  ERROR  * 

1 / 1 0 > » 1 6 F PRO GRAM  STOPPED  ,10!/)) 

504  FORMAT (10!/) .10X.40HI NPUT  ERROR  IN  EXCEPTIONS  TO  SHAFT  DATA 
1 / 1 0 X , 16  FPROGRA  M STOPPED  ,10!/)) 

505  FORMAT ( 1 0!/) , 10X, 37HINPUT  ERROR  IN  UNIT  TYPE  DESIGNATION  / 

1 1 OX , 16FPR0GRAM  STOPPED  ,10!/)) 

5C6  FORMAT (10!/), 10X, 37HINPUT  ERROR  NO.  OF  FLOORS  EXCEEDS  NH  / 

1 10X  « 1 6H FRO GRAM  STOPPED  .10!/)) 

5 C 7 FORMAT! 10!/) . 10X .25HINPUT  ERROR  IF1  .GT.  IF2  / 

1 1 OX  , 1 6HFR0GRAM  STOPPED  ,10!/)) 

510  FORMAT (10!/) « 10X, 36HI NPUT  EXCEEDS  DIMENSION  PARAMETER  , A3/ 

+ 1 0 X , 1 6F  PROGRA  M STOPPED  ,10!/)) 

530  FORM  AT ! 1 OX ,3A6) 

535  FORMAT!//  10X.26HFL0W  COEFFICIENTS  AS  READ  ) 

540  FORM  AT ( 1 OX, 15HBUILDING  DATA  //34X , 1 1 HTE MPERATURE  .4X.5HFIXED, 

1 12X  ,2! 1 IX , 4HFL0 W ) / 1 OX , 1 1H CO MPARTMENT  , 4X , 5HFL00R , 6X, 7HPR0F ILE . 

2 6X.4FFLQW,  5X  « 1 3HCCNNECT  ION  TO  » 4X  , 11  FiCOEFFIC  I ENT  ,4X, 

3 8 F AREA  ) 

541  F0RMAT(/4X.3I12.F14.i ) 

5 42  FORMAT  !/4X. 3112. F 14.1 .4X.5HP0INT ,I7»F 11 .2 ,F1 5.4) 

5 43  FORMAT (58X  ,5HPC INT , I7.F11.2.F15.4) 

544  FORMAT (58X,9HOUTS IDE  • I 3 . F 1 1 . 2 ,F 1 5. 4 ) 

END 
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SUEROUTINE  CORR 


£NBS*PLI  E$  .SEC IN  A. CORR 
SUERCUTINE  CORR 

THIS  ROUTINE  CALCULATES  ADJUSTED  FLOW  COEFFICIENTS 
(C1.C2.CC1 . CC2 ) 


PARAMETER  (MM  = 14  0 ,MS=e . MC=9 , MFG=2 , MT P=2 , MFL=25 . MB=50 ) 

CCMMCN  /CORR/C1 (MM , MC ) . C2 (MM , MC ) .CO  1 ( MM , MPO ) . C02(MM,MP0 ) 

COMMON  NT.  P( MM  ) ,C( MM.MC  ) .NC(MM)  , JC( MM, MC) . ITS(MS) , 

1 FC( MM  ,MC)  ,PZ( MM, MC) ,PO( MM.MPO ) ,CO( MM, MPO) ,F(  MM) , PFO(MFL.MPO ) , 

2 FF(MM),FO( MM , MFO  > ,CS( MS) . PS ( MFL ) .NS1 (MS) ,NS2(MS) . 

3 FSS (MS  ) .N.NS  ,NPO,  ICONV ,E, I BUG .AI (MM, MC) . AO( MM, MPO) , TITSH(MS .5)  , 

4 NH, M MFL) , IFLOOR ( MM ) , T ( MTP , MFL ) .NFSl ( MS) ,NFS2( MS) , IT( MB) ,NTP 

5 ,NCC( MM ) . JOC( MM, MFC) , TOUT 

COUELE  PRECISION  P.FO.PS 
CO  12  1=1. NT 

CORRECT  C 

PATMCS= 101325. 

EB=  1 CO 0.*SQRT(2.*PATMO S/287. )/ 1.2 
NN=NC(  I) 

IF  ( I . CT.  N ) GO  TO  1 
IP=IK  I) 

CO  TC  4 
CO  2 I S= 1 , NS 

IF ( I .LE.  NS2 (IS)  .ANC.  I . GE . NS1(IS))G0  TO  3 
COM  INUE 
WRITE(6»700) 

STOP 

ip=nsas> 

IFF=  IF  LOGR  ( I) 

T1 =T ( IF, IFF) 

IF (Nh  .EQ.  0 ) GO  TO  10 
CO  9 J = 1 ,NN 
JJ  = JC ( 1, J) 

Cl  (I  , J )=BB*C(  I, J) *AI ( I , J ) /SORT ( T1 ) 

IF ( J J .GT.  N) GO  TO  5 
IP=IT ( JJ ) 

GO  TC  8 
DO  6 I S=  1 , NS 

I F ( J J .LE.  NS2 ( IS ) .ANC.  JJ  .GE.  NS  1 ( I S ) ) GO  TO  7 
CCNT  INUE 
WRITE  (6,700) 

STOF 

IP=ITS( IS) 

IFF  = IFLOOR( JJ) 

T2=T  ( IP,  IFF) 

C2 ( I , J )=BB*C(  I . J) *AI ( I . J)/SQRT(T2) 

CONT  INUE 

CORRECT  CO 

10  NN  C = M CO  ( I) 

IF ( NhC  .EO.  0 ) GC  TC  12 
CO  11  J=  1 , NNC 

C0  1 ( I,  J)=BB*CC( I, J )*AG( I , J)/SGRT(T1 ) 
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CO  2 ( 1. J)=6B*C0< It J )*AC< ItJ)/SQRT(TOUT) 

11  CONTINUE 

12  CONTINUE 
RETURN 

700  FORM  AT (/// 1 OX  » 36HPR0GRAM  STOPPED  IN  SUBROUTINE  CORR 

END 
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SNBS*PLIES  .SHCfc  A.  INIT 
SUERCUTINE  INIT 
C 
C 

C THIS  ROUTINE  INITIALIZES  Tl-E  PRESSURE  ARRAY 

C 

PARAMETER  <MM= 140 , MS  = 8 , MC=9, MPG=2 , MTP=2 ,MFL=25 , MB=50) 

PARAMETER  (MEP=MB+1) 

COMMON  /C0RR/C1 (MM  ,MC)  . C2 ( MM , M C ) * CO  1 ( MM , M PO ) , C02( MM.MPO ) 

COMMCN  NT.  P(MM  ) ,C(MM,MC)  ,NC(MM)  , JC ( M M , MC ) . ITS(MS)  . 

1 FC( KM . MC) .PZ ( MV, MC) .PC( KM, MPC ) ,CO( MM , MPO ) ,F( MM ) , PFO( M^L, MPO ) , 

2 FF ( MM  ) ,FC( MM.MPO ) ,CS( MS) »PS(MFL) ,NS1 ( MS)  , NS2 ( MS)  , 

3 FSS  (MS)  ,N,NS  , NFO , ICCNV.E, I EU G . AI ( MM , MC ) .AO(MM.MPO)  ,TITSH(MS,5)» 

4 NH , K MFL ) . IFLOCR ( MM)  ,T< MTP, MFL) .NFS1 ( MS)  ,NFS2( MS ) . IT(MB  ),NTP 

5 ,NCC ( MM  ) , JOC( MM, MPC) , TOUT 
COUELE  PRECISION  P.PO.PS 
CIMENS1CN  SC  ( MS  ) , SCO  (MS) 

COMMCN  /MAT/A( ME. MEP) ,XX ( MB) ,NNN 
COUELE  PRECISION  A, XX 

NNN  = N 
C 

CALCULATE  AVERAGE  OUTSIDE  PRESSURE 
SUM=C  . 

CO  1C  J=l,NPO 
CC  1C  1=1, NH 
SUM  = SU  M-f  PFQ( I , J) 

FA=SLM/( NPO*NH) 


C 

C Tt-E  CO  LOOP  TO  STATEMENT  30  ESTIMATES 

C SHAFT  PRESSURES 

C 

DO  30  IS  = 1 , NS 
C 

C CALCULATE  SHAFT  PRESSURE  C I F FERENCE  , DP 

C 

SUM  = C . 

SUMN  = 0 . 

N 1 =N  S 1 ( I S ) 

N2=NS2 ( IS) 

CO  16  I=N1 , N2 
SUM=SUM+FF( I) 

NN  =NC ( I) 

IF  ( N N .EC.  0 . ) GO  TO  16 
DO  IE  J=  1 , NN 
SUMN  = SUMN+CK  I , J) 

15  CONTINUE 

SC ( IS  )=SUMN 

16  NN  C = NCC(  I) 

IF (NNO  .EQ.  0 ) GO  TO  18 
DO  17  J=  1 ,NNO 
SUM,N  = SUMN  + C01  ( I , J ) 

17  CONTINUE 

SCO ( IS  ) = SUMN—  SC ( I S ) 

18  CONTINUE 
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SUBROUTINE  INIT 


CP 2=  SUM/SUMN 
SI GN  = 1 . 

IF  C CF2  .LT.  0.)SIGN=-1. 

CP  = S 1GN*CSIGN*DP2)**2 
C 

C CALCULATE  AVERAGE  TEMP  CF  SHAFT 

C 

SUM=C. 

ip=itscis) 

DO  20  I = N 1 • N2 
IFF  = IF LOCR(  I ) 

20  SUM=SUM+TC IP,  IFF) 

TA  = SLM/C  N2-NI  + 1 ) 

C 

C ESTIMATE  PRESSURE  AT  BCTTCM  OF  SHAFT  . 

C 

HH=0.5*CHCNH)-HC1 ) )+HCl  ) 

NF 1 = N FS 1 ( IS) 

PB0T=PA+DP+3462.*C HH-F(NF1 ) )/TA 
C 

C ESTIMATE  OTHER  SHAFT  PRESSURES 

C 

PCM  ) = PBCT 
NM  =N  2 - 1 
CO  24  I=N1,NM 
IP1=1+1 

24  PC  IP  1 ) =P( I )-PZ( I. 1 ) 

3 0 CONTINUE 

C 

C ENC  CF  SHAFT  PRESSURE  ESTIMATES 

C 

C SET  UP  MATRIX  FOR  BUILDING  COMPARTMENTS 

C 

NP  1 = l 
CO  50  1=1,  N 
NN  = NCC  I) 

SUMI  1 = 0. 

SUMNF=0. 

IFCNN  .EC.  0 . ) GC  TO  42 
CO  4 C JJ=1,NN 
J= JCC  I ,J J) 

IF  C J .CT.  N ) GO  TO  34 
AC  I ,.)=C1C I.JJ) 

SUMI  1=SUMII-C1C  I.JJ) 

SUMNF=SUMNP— Cl C I, JJ)*PZC I.JJ) 

GO  TC  40 

34  SUMI  1=SUM1I-C1C  I.  JJ> 

SU  MNF=  SUMNP— C 1 C I*  J J ) *P  C J ) 

40  CONTINUE 

42  NNC=NCCC I ) 

IF CNNO  .EQ.  0 ) GO  TO  46 

CO  45  K= 1 , NNO 

SUMI I = SUM I I— CO  1 C I , K) 

45  SU  MN  F = SUMNP— CO 1CI  »K)*FOCI»K) 

46  AC  I,  I ) = SUM 1 1 

AC  I * NP  1 ) = SUMNP— FF  C I) 

50  CONTINUE 

C 


PBOT 
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c*  r.  n 


SUBROUTINE  INIT 


C WF  ITE  MATRIX 

C 

IF  < I EU 0 .EQ.  0 ) GO  TO  e4 
V»RITE(€.802) 

CO  £2  1=1.  N 

£2  WRITE(£,803)(A( I, J)  .J=l ,NP1> 

C 

c 

C CALL  ROUTINE  TO  SOLVE  FOR  INITIAL  BUILDING  PRESSURES 

C 

84  CALL  SJMEQ 

C 

C OUTPUT  INITIAL  PRESSURES 

C 

IF ( I EU  G .EQ.  0 ) GO  TO  85 
WRITE(  8,800) 

WRITE  (6, 801  )(  I,XX<  I).  1 = 1 ,N> 

NN=NS1 ( 1 ) 

WR  ITE  (6. 801 )(  I.P(  I ) « I=NN . NT ) 


ASSIGN  EUILDING  PRESSURES 
C 


89 

SO 

DO  SC  1=1, N 

pc i >=xxc  i) 

RETURN 

800 

801 

FORMAT C///8C6X , 1HI  .4X.3HP 
FORM  AT (8(17. F7.1)  ) 

)/) 

802 
6 C3 

F OR MAT( /// 10X, 20 H MATRIX 
FC  RM  AT  ( 1 OX  . 1 1 F 1 1 . 1 ) 

ENC 

COEFFICIENTS 

/ ) 

a ♦"  DG  • 

P 

SUBROUTINE  BLDGP.L.l 
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A A A A A A A A A A M A A A A A A 


SUBROUTINE  BLDGP 


ENBS*PLI E$ .SHOW  A. BLDGP 
SUERCUTINE  BLDGP 
C 

c 

C THIS  ROUTINE  CALCULATES  STEADY  STATE  PRESSURES 

C FOR  BUILDING  COMPARTMENTS 

C 

c 

c 

PARAMETER  (MM=140  , MS= e , MC=9 , MPC  = 2 . MTP=2 , MFL=25 , MB=50 ) 

COMMON  NT,  P(MM ) ,C( MM ,MC ) ,NC (MM) . JCC MM, MC > , ITS( MS > , 

1 FC( MM  ,MC) ,PZ ( MM,  MC) , PO( MM, MPO ) , CO ( MM ,MPO ) ,F( MM ) , PFO ( MFL , MPO ) * 

2 FF  (MM  ) , FO  ( MM  , MPO  ) ,CS(  MS)  ,PS(  MFL)  ,NS1  ( MS)  ,NS2  (MS)  , 

3 FSS(MS) ,N,NS • NFO , I CON V , E , I BUG , AI ( MM , MC ) , AO( M M, MPO ) , TI T SH (M S ,5 ) , 

4 Nh.KMFL),  IFLOOR  ( MM  ) , T ( MTP  , M FL  ) .NFS!  ( MS  ) ,NFS2(  MS ) . IT < MB ) , NT P 

5 , NCC( MM ) , JOC( MM, MFC)  .TOUT 
CCUELE  PRECISION  P.FO.PS.PI 
IF ( I E U G .GT.  0 ) WR  I TE(  6,806 ) 

IT  M = 100 

ICON V=0 
CO  IS  1=1,  N 
C 

C CALCULATE  NET  FLOW  .FI.  INTO  POINT  I 

FI=PFLCW ( I ,P( I ) ) 

CHECK  MAGNITUDE  OF  FI 
IF  (AES(F  I)  .LT.  E ) GO  TO  IS 
ICCNV=ICCNV+1 

SET  UP  PARAMETERS  FOR  ITERATION 
CP  = 1 .0 
IPHA  SE  = 1 
CP  I = C. 

EE  = 0 .2*ABS(FI ) 

IF  ( EE  .LT.  E ) EE=E 
SI GN  = 1 

IF ( F 1 .LT.  0.  >SIGN=-1 
IK  = 0 

IF ( I EU  G .GT.  0 ) MR  I TE (6.602) 

ITERATION  TO  RECUICE  MAGNITUDE  OF  FN 
IK  = I K+  1 

NEW  ESTIMATE  OF  PRESSURE  .PI.  AT  POINT  I 
PI  = F ( I ) + S I GN*QP 

CALCULATE  NET  FLOW  , F N , INTO  POINT  I USING  PI 
FN  = F FLCW ( I , P I ) 

IF ( IEUG.GT.O) WRITE! 6.8  04)1, IK , F I , FN , FP , DP  I , DP , DPP , P I , I PHASE 

CHECK  MAGNITUDE  OF  FN 
IF(AES(FN)  .LT.  EE ) GO  TO  10 

CHECK  NUMBER  OF  ITERATIONS 
IF(  IK  .GT.  I T M ) GO  TO  25 

CHECK  PHASE 


47 


<u  r>  r> 


SUBROUTINE  ELDGP 


IF  ( IFHASE  .EQ.  2)  GO  TO  6 
C 

C CHECK  FOR  TRANSITION  FROM  PHASE  1 TO  PHASE  2 

IF ( F I *FN  .LT.  0.)GC  TO  A 
C 

C PHASE  1 

CP  I = CP 
CP  = 5 .0  *DP 
F I =F  H 
C-0  TC  2 
C 

C PHASE  2 

4 IPHASE=2 
GO  TC  S 

t IF ( F I*  FN  .GT.  0.)GC  TO  8 

C 

C NEW  CP  EETWEEN  DPI  AND  DP 

5 CPP=CP 
FP  = FH 

CP  = DF !♦<  CPP-DP I )*F  I/(  FI-FN ) 

GO  TC  2 
C 

C NEW  CP  EETWEEN  CP  AND  CFF 

€ FI  = FN 

CP  I = CP 

CF=DFI4< CPP-DP I )*FN/(FN-FP) 

GO  TC  2 
10  P(  I ) =P  I 

15  CONTINUE 

C 

RETURN 

25  WR  I T E ( 6 t 80  0 ) 

STOF 

C 

FORMAT  STATEMENTS 

00  FORMAT  (///10X, 20(  1 H* >/// 1 OX, 22HEXCESSI VE  ITERATIONS  / 

+ 10X.8HIN  ELDGP  ///10X»20<1H*) //////) 

£02  FORMAT (//l IX, 1HI , 2X , 2H I T . 12 X, 2HF I , 1 3X « 2HFN*  13X, 2HFP, 12X  , 3HDPI, 

+ 13X, 2H CP.  1 2X. 3H DPP* 13 X, 2 HP I ,3X,5HPHASE  /) 

£04  FORMAT  (8X,2I4,3E1 5. 4, 4F 15. 6. 15) 

£06  FORM  AT ( /// 1 OX, 6HBLDGF  ) 

END 


aFDG.P 


SUBROUTINE  SHAFTP.L.l 


n n n n n n n n m n n n a n a a a a a a a a 


SUBROUTINE  SHAFT® 


aNES*PLI  ES  . SFO  I*  A.SHAFTF 
SUERCUTINE  SFAFTP 


C 


THIS  F 0 U T I NE  CALCULATES  STEADY  STATE  PRESSURES 
FCR  SHAFTS 


PARAMETER  ( M M= 14 0 . MS=8 • MC  = 9 » MFC  = 2 , MTP=2 . MFL=25, MB=50 ) 

COMMON  NT.  F<  MM  ) ,C(  MM  ,MC  ) ,NC(MM)  ♦ JC<MM,  MC  > . ITS(MS)  . 

1 FC(MM,MC),PZ(MM,MC),PO(MM,MPC) ,CO  < MM  , MPO  ) , F ( MM  ) , PFO C MFL . MPO ) , 

2 FF(MM)»FO( MM . MPO ) .CS( MS) »PS(MFL) .NS1 ( MS) .NS2(MS) . 

3 FSSCMS)  .N.NS  . NPO . ICONV.E, IBUG.AI (MM.MC) . AO( MM*  MPO )»TITSH(MS,5)  . 

4 Nt-.KNFL),  IFLCCR  ( MM  ) ,T(  MTF.  MFL)  .NFS1  ( MS)  .NFS2(  MS  ) . I T ( M B ) .N  TP 

5 ,NCC(MM), JOCCMM, MPO). TOUT 
OOUELE  PRECISION  P.PO.PS.PI 
IF ( I EU  C .GT.  0 )WR ITE( 6,806) 

ITM=  100 
ICGN  V = 0 
DO  IE  1=1, NS 

CALCULATE  NET  FLOW  ,FI,  INTO  POINT  I 
N1=NS  1 ( I ) 

FI=SFLCW ( I ,P( N1 ) ) 

CHECK  MAGNITUDE  OF  FI 
IF(AES(FI)  .LT.  E ) GO  TO  15 
ICON  V = ICCNV  + 1 

SET  UP  PARAMETERS  FOR  ITERATION 
CP  = 1 .0 
IPHA  SE  = 1 
DP  I=C. 

EE  =0  «2*ABS( FI ) 

IF ( EE  .LT.  E) EE=E 
S I GN  = 1 

IF  (F  I .LT.  0.  )SIGN=— 1 
IK  =0 

IF ( I E U G .GT.  0 )WR ITE( 6 ,802 ) 

ITERATION  TO  REDUICE  MAGNITUDE  OF  FN 
IK  = IK+  1 

NEW  ESTIMATE  OF  PRESSURE  .PI,  AT  BOTTOM  OF  SHAFT  I 
PI  = F ( N 1 )+SIGN*DP 

CALCULATE  NET  FLOW  , F N , INTO  SHAFT  I USING  PI 
FN=  S FLOW ( I ,P I ) 

IF(IEUG.GT.0)WRITE(6,804)I.IK,FI,FN.FP,DPI , DP , DPP, PI , I PHASE 

CHECK  MAGNITUDE  OF  FN 
IF (AES (FN  ) .LT.  EE ) GO  TO  10 

CHECK  NUMBER  OF  ITERATIONS 
I F ( I K .GT.  IT  M ) GO  TO  25 
C 

C CHECK  FHASE 
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r\  n in  r>  n ft 


SUBROUTINE  SHAFTP 


IF  ( I F F AS  E .EQ.  2)G0  TO  6 
C 

C CFECK  FOR  TRANSITION  FROM  PHASE  I TO  PHASE  2 

IF(FI*FN  .LT.  0 • ) G 0 TO  4 
C 

C FHASE  1 

CP I=CP 
CP  = 5 .0*0P 
FI  =Fh 
GO  TC  2 
C 

C FHASE  2 

4 IPHA  SE  =2 

GO  TC  9 

IF ( F 1*  FN  .GT.  0 • ) GO  TO  8 

NEW  CP  BETWEEN  CPI  AND  DP 
CP P=CP 
FP  = FN 

DP=DFI+( CPP-DPI )*FI/(FI-FN) 

GO  TC  2 

NEW  CP  E ET WEEN  DP  AND  DPP 
E F I =F  N 

CP I=CP 

DP=CFI+( DPP— DPI )*FN/(FN-FP> 

GO  TC  2 

10  N2  = N £2(1) 

CO  11  IF=N1,N2 
II = IF+  1-N1 

11  F(  IF  1 = PS (II) 

15  CONTINUE 

C 

RETURN 

25  WRITE(C.800  ) 

STOP 

C 

C FORMAT  STATEMENTS 

C 

€00  FORMAT (///10X »20( 1 H* > /// 10X , 22HEXCESS I VE  ITERATIONS  / 

+ 1 OX .9  H IN  SHAFTP  /// 1 OX , 20 ( 1 H* )////// ) 

€02  FORMAT {//l IX.  I H I , 2 X . 2H I T • 12X . 2HF I . 1 3X • 2HFN . 1 3X . 2HFP . 1 2X . 3HD  PI. 

+13X.2HCP  . 1 2X . 3FDPP  .13X.2HPI .3X.5HPHASE  /» 

€04  F0RMAT(8X,2I4,3E1 5.4.4F15.C .15) 

€06  FORMAT ( ///10X, CHSHAFTP) 

END 


£ H DG . P 


SUBROUTINE  PZAD 


L.l 


SUBROUTINE  PZ AD 


£NBS*PLI  E$  .SHOW  A.PZAD 
SUERCUTINE  PZAD 
C 

C THIS  ROUTINE  CORRECTS  PZ  TERMS  FOR  PRESSURE 

C 

PARAMETER  ( M M= 1 40  . MS=8 . MC  = 9 , MFC=2 , MTP=2 , MFL=25 , MB=50 ) 

CQMMCN  NT,  P ( MM ) , C { MM , MC ) , NC ( MM ) , JC ( MM , MC ) , ITS( MS  ) , 

1 FC( MM  ,MC>  ,PZ ( MM, MC) ,P0( MM ,MPC) ,CO( MM , MPO ) , F(  MM  ) , PFO ( MFL  » MPO ) , 

2 F F ( M M ) , FO( MM , MPO ) ,CS( MS), PS (MFL) ,NSl ( MS)  ,NS2(MS)  , 

3 FSS  (MS  ) ,N,NS  ,NFO , ICONV, E. I EUG , AI ( MM , MC ) , AO( MM, MPO) ,TITSH(MS,5). 

4 NH, F( MFL) , IF LCCR ( MM) , T ( MTP , MFL) ,NFSI  ( MS) ,NFS2(MS), IT(MB)»NTP 

5 ,NCC (MM ) , JOC( MM, MPO) , TOUT 
COMMCN  /PZ Z/  PGZ 

COUELE  PRECISION  P,F0,PS 
IF ( I EUG  .GT.  — 2 ) G O TO  1 
WRITE (€.800) 

DO  2 1 = 1 ,N 
NN  = NC(  I ) 

IF  ( N N .EQ.  0 ) GO  TO  2 

WRITE  (€.801  ) ( I.J.PZd,  J),J=1.NN) 

2 CONTINUE 

NP 1 = N + 1 

WR  ITE ( €, 802 )( IL,PZ( IL, 1 ) ,IL  = NP 1 , NT) 

1 CO  1C  1=1, N 

NN  = NC(  I) 

I F ( N N .EC.  0 ) GO  TO  10 
IA= IT ( 1) 

IF  I = IF LOCR( I ) 

DO  8 JJ=1,NN 
J=  JC  ( I , J J) 

IF  J=  IFLOOR  ( J ) 

I F ( IFI  .EQ.  I F J ) GO  TO  8 
IB= IT ( J ) 

TE  MP A = 0 .5*(T(IA,IFI)  + T(  I B , I F J ) ) 

FAVE  = 0 .5*  (P  (I  )+P  ( J ) ) +PGZ 

PZ ( I «JJ)=(0. 03416  *PA VE/TEMPA )* (H( IF  J ) — H ( IFI ) ) 

€ CONTINUE 

10  CONTINUE 

CO  20  IS=  1 , NS 
N l =N  £ 1 ( I S ) 

N2=NS2(IS)-1 
ITT=  1TS(  IS) 

CO  IS  I=N 1 ,N2 
IF  I = IFLOQR(  I ) 

IF  J = IF  1+  1 

TEMPA=O.S*(T(ITT,IFI)+T( ITT. IFJ) ) 

J=  1+  1 

PA=  0 • 5 * ( P( I ) +P ( J ) ) +PGZ 

15  PZ ( I , 1 )=(0.034l€*PA/TEMPA)*(H( IFJ)- H( IFI) ) 

20  CONTINUE 

RETURN 

E00  F0RMAT(/10X ,1 OHINI TIAL  FZ  /) 

601  FORMAT ( 10X ,3HPZ( , 12, 1H, 12. 4H ) = ,F12.4) 

£02  FORM  A T ( 1 OX , 3HPZ (♦  I2.6H.1  ) • ,F12.4) 

£03  F0RMAT(/10X  ,HHAD  JUSTED  PZ  /) 

END 
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SUBROUT  INE  OUT 


aNES*PLIE$ .SHOW  A. OUT 
SUERCLTINE  OUT 
C 

c 

C THIS  ROUTINE  OUTPUTS  FLOWS  AND  DIFFERENTIAL  PRESSURES 

C FOR  ALL  SHAFTS  ANO  BUILDING  COMPARTMENTS 

C 
C 

PARAMETER  ( MM= 14 0 , MS=8 . MC=9 . MPC=2. MTP=2 . MFL=25 . MB=50 ) 

CO  M VC  N /CORR/C1 (MM  ,MC  ) ,C2 (MM , MC ) ,COi ( MM, MPO) » C02( MM, MPO ) 

COMMCN  / IO/T ITLE ( 18)  , IOUT , I UN I T , NCOMP ( MFL ) ,SNCOMP(MFL) 

COMMON  NT,  P( MM  ) ,C(  MM , MC)  ,NC( MM) , JC( MM.MC) , ITS(MS ), 

1 FC ( MM ,MC) ,PZ( MM, MC) , PO ( MM, MPO ) ,CO ( MM , MPO ) , F(  MM) , PFO ( MFL , MPO ) , 

2 FF ( MM)  ,FO( MM.MFO)  ,CS<  MS) , FS( MFL ) *NS1  ( MS) ,NS2 (MS ) , 

3 FSS  (MS  ) ,N,NS  ,NFO , ICONV ,E , I BUG  . AI (MM.MC), AO (MM, MPO) , TITSH(MS .5 ) , 

4 NH,  H(  MFL ) , I FLOOR ( MM) , T( MTP .MFL ) ,NFS1 ( MS)  , NFS2( MS) , IT (ME) ,NTP 

5 , NCC ( MM  ), JOC( MM,  MFC) .TOUT 
CCUELE  PRECISION  P.FC.FS 
INTECER  COM 

C 

IUN  IT  = 1 FOR  SI  UNITS 
I U N I 7 = 2 FOR  ENG  UNITS 
WHEN  I UN IT  = 2 GO  TO  100 
IF ( I UN  IT  .EQ.  2 ) GO  TO  100 

BUILDING  COMPARTMENT  OUTPUT 


1 = 0 

IL  = 0 

WRITE(e,800)(TITLE(  I)  .1  = 1.18) 

CO  2 C 1FF= 1 . NH 
NNN=NCCMF( IFF) 

IF (NNN  .EQ.  0 ) GO  TO  30 
DO  29  IC  =1 .NNN 
1=  14  1 
NN  =NC ( I) 

NNC=NCC ( I) 

IL  = I L 4 HN  4NN04  2 

IF  ( IL  .LT.  £ 1 ) GO  TO  2 

taR  I T E ( €.80  0 )( TI TLE ( I ) ,1=1 , 18 ) 

IL  = NN4NN042 

IF ( N N .GT.  0 ) GO  TO  3 

WR  ITE ( 8, 801 ) IFF, IC,P(I ) , I T ( I ) ,FF(  I ) 

GO  TC  21 
CO  2 C J=  1 , NN 
JJ= JC ( I . J) 

CP=  F (JJ)-P(I)4PZ(  I.J) 

CC  = C 2 ( I.J) 

IF ( D F .LT.  0. )CC=C1 ( I , J) 

IF  ( J J .LE.  N ) GO  TO  10 
CO  5 I S=  1 . NS 

I F ( J J .GE.  NS  1 ( IS  ) .AND.  JJ  .LE.  NS2(IS))G0  TO  6 
CONTINUE 

IF  ( J .C-T.  1 ) GO  TO  7 

WR ITE(C, e02) IFF, IC .P ( I ) . IT ( I ) . FF ( I ) , ( TI TSH( IS .K ) ,K=1 ,5) 
4 , DF ,CC. AI  ( 1,1)  ,FC  ( 1,1) 

GO  TC  20 
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SUBROUTINE  CUT 


7 WR  IT  E ( £ . 603)(TITSH( I S . K ) « K= 1 ♦ S ) .DP  . CC • A I ( I.J)  .FC(I.J) 

C-0  TC  20 

10  IF  J = IF  LOOK  ( J J ) 

CO*=JJ-SNCOMP( I F J ) 

IF  ( J .GT.  1 )G0  TO  12 

WRITE!  £.804)  IFF,IC.P( I).IT(I),FF(I) , I F J , COM. DP. CC * A I ( I . 1 ) .FC(  I • l ) 
GO  TC  20 

12  WR  ITE(6.805)IFJ,CGM»DP.  CC  * A I ( I ,J),FC(I. J) 

20  CONTINUE 

21  IF ( NNG  .EQ.  0 ) GO  TO  29 

CO  2 2 J-l.NNO 

JJ  = JCC (I.J) 

CP=PC( I * J )-P( I ) 

CC=CC2 (I.J) 

IF ( 0 F .LT.  0 . ) CC=C01 (I.J) 

23  WRITE! £.806) JJ.DP .CC • AO (I.J) *F0( I.J) 

29  WR  ITE  (£.  807  )F(  I ) 

30  CONTINUE 

WR IT  E ( £ • 900 ) 

C 

C SFAFT  OUTPUT 

C 

CO  6 C IS  = 1 . NS 
N1=NS1 ( IS) 

N2=NS2 (IS) 

22  WRITE(£.ei4)(TITLE(I)  .1  = 1.18) 

WR IT  E ( £ . 80  8)(TITSH(IS*K).K  = 1,E)  . ITS( I S ) *CS( IS ) 

CC  SC  I=N1 • N2 
NN  = NC(  1) 

IF( NN  .GT.  0 ) GO  TO  35 

WR  ITE(£. 609)1 FLOOR (I).P(I)  .FFC I) 

GO  TC  41 

25  CO  4 C J=  1 . NN 

JJ=JC(  I.J) 

CP=P  ( J J )— P ( I ) 

CC  =C  2 ( I.J) 

IF  ( OF  .LT.  0.  )CC=C1(I.J) 

IF  J=IFLOCR(  JJ) 

COI»=  JJ-SNCOMP(  IFJ  ) 

IF  ( J .C-T.  1 ) GO  TO  36 

WR  ITE(6. 810)  I FLOOR  (I).P(I).FF(I)  , I F J,  COM  , DP  ,CC  * AKI,  l ) .FC(I.l) 

GO  TC  40 

2£  WRITE  (£.81  1 )IFJ«CCM»DP»CC»  AI  (I.J).  FC  (I.J) 

40  CONTINUE 

41  NNC  = NCC ( I ) 

IF (NNO  .EQ.  0 ) GO  TO  50 
DO  4£  J= 1 . NNO 
JJ  = JCC  (I.J) 

CP  = FC  ( I. J )-P(  I ) 

CC  =CC2 ( I .J) 

IF  ( CF  .LT.  0 • ) CC=CO  1 ( I . J ) 

46  WRITE(6.812)JJ.DP,CC.A0(I,J),FG(I,J) 

£0  CONTINUE 

WRITE(£ ,813)FSS( IS) 

WR  ITE (£.900 ) 

£0  CONTINUE 

GO  TC  165 
C 
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SUBROUTINE  OUT 


C EUILCING  CATA  OUTPUT  FOR  IUNIT  = 2 

C 

100  1=0 
IL  = 0 

WR  ITEtC.eOO )( TITLE! I)  .1=1  . 18) 

CO  120  IFF=1,NH 
NNN=NCCMF { IFF  ) 

IF ( NNN  .EQ.  0 ) GO  TO  130 
CO  129  IC=1 »NNN 
1=1+  1 

FF  1=  F ( 1 )/0 . 47  19 
PI II=P(I )/248.8 
FF  F = FF ( I )/0.4719 
NN=NC ( I ) 

NNC=  N CO  C I) 

IL  = I L + NN+NNO+2 

IF  ( IL  .LT  . 51 ) GO  TO  102 

WR ITE( €.80  0 )( TI TLE(  I ) .1  = 1 , 18  ) 

IL  = N N +NN0+2 

102  I F ( N N .GT.  0 ) GO  TO  103 

WR  ITE( €, 601 ) IFF, I C. PI 1 1 , IT< I > »FFF 
GO  TC  121 

103  CO  120  J = 1 » NN 
FCCC=FC(  I, J )/0 .47 1 9 
JJ  = JC  ( I , J) 

CP= ( F ( J J )-P(I )+PZ ( I,J»  )/24  8.8 
AA 1=41 ( I » J ) /O • 092  9 
CC  = C 2 ( 1. J) 

IF ( D F .LT.  0. )CC=C1( I.J) 

CC=CC*33  .43 

IF ( J „ .LE.  N ) GO  TO  110 
CO  1C5  I£=l .NS 

IF ( J J .GE.  NS  1 ( IS ) .AND.  JJ  .LE.  NS2(IS))GO  TO  106 

105  CONTINUE 

106  IF  ( J .GT.  1 ) GO  TO  107 

WR ITE<  6.60  2) IFF.IC.PIII . ITC I I .FFF  »(TITSH(  IS  «K  ) , K= 1,5) 

+ .CP.CC.AAI, FCCC 
GO  TC  120 

1 C 7 WR  IT  E ( € , 602  )(TITSH(  IS.K  ) ,K=1 .5)  .DP  ,CC  ,AAI  .FCCC 

GO  TC  120 

110  IF J= 1 FLO OR ( JJ  ) 

CO  ^ = J J - S NCC  MP  ( IF  J ) 

IF  C J . GT  . 1 ) GO  TO  112 

WR  ITE(£,604)IFF  .IC.PIII , IT  (I ) » FFF  , I F J , COM, DP , CC . AAI . FCCC 
GO  TC  120 

112  WRITE! 6, €05) I F J ,CC  * » DP ,CC ,AAI .FCCC 

120  CONTINUE 

121  IF ( NNO  .EQ.  0 ) GC  TO  129 
CO  123  J = 1 » N NO 

FOC  = F 0 ( I , J )/0 • 471 9 
JJ  = JCC ( I • J ) 

DP= (FC  (I.J  )-P(  I ))/248.8 
AAC=  AC ( I , JJ/0.0929 
CC  = CC  2 ( I.J) 

IF ( D F .LT.  0 • ) CC  = C0 1 (I.J) 

CC=CC*22 .43 

123  WRITE(€.€06) JJ.CP  »CC , A AC , FOO 

129  WR  ITE ( €, 807)FF I 
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SUBROUTINE  OUT 


130  CONTINUE 

*R  1TE(£.901) 

C 

C SI-AFT  OUTPUT  FOR  IUNIT  = 2 

C 

DO  ISO  IS=1,NS 
CSS=CS< IS)/0. 02992 
FF  1=  F S S ( IS)/0.4719 
N1=NS1  (IS) 

N2=N£2 ( IS) 

132  IaR  I T E < £,814)( TITLE! I ) «I=1 . 18) 

IaR  IT  E ( £ , 80  8 ) ( T ITSH  < I S . K ) , K=1  , 5 ) * IT  S<  I S)  .C  SS 

DO  ISO  I=N1,N2 

FFF=FF ( I )/0.4719 

FI  I I=P(I  )/248 . 8 

NN=NC(  1) 

IF ( N N .GT.  0 ) GO  TO  135 

IaR  ITE(£.€09)  I FLOOR  ( I ) * PI  II  »FFF 

GO  TC  141 

135  DO  1 A 0 J=1,NN 
FCCC=FC(  I • J )/0 .4719 
JJ=JC ( It  J) 

DP  = ( F ( J J ) — P C I ) )/24  8.8 
A A 1=  A I ( I , J )/0  *0929 
CC=C  2 ( I» J) 

IF ( C F .LT.  0 • )CC=C 1(1 t J ) 

CC=CC*23.43 
IF J=  IF  L0CR( JJ) 

CCM= JJ-SNCOMP!  IF  J ) 

IF ( J .GT.  1 ) GO  TO  136 

WRITE(€,610) I FLOOR (I).PIIItFFF  , IF J , COM . OP ,CC , AAI .FCCC 
GO  TC  140 

136  <aRITE(£,€11)IFJ,C0M,DP,CC,  AAI, FCCC 

140  CONTINUE 

141  NNC=NCO(  I) 

IF ( NNO  .EQ.  0 ) GO  TO  150 
CO  146  J= 1 , NNO 
FOC=FO(I , J ) /0 . 471 9 
JJ=JCC ( I , J ) 

CP  = (FO(I  ,J)-P( I }) / 24 8.8 
AAG  = AO ( I , J )/0.0929 
CC=CC2 ( I ,J) 

IF ( CP  .LT.  0. )CC=C01( I . J ) 

CC  =C  C*  23  .4  3 

146  \aR  ITE(£,£12)JJ,  DP  , CC  , A AO  , FOO 

150  CONTINUE 

IaR  ITE  ( £.  813  )FF  I 
IaR  I T E ( € , 90  1 ) 

ICO  CONTINUE 
C 

C SUMMARY  OUTPUT 

C USER  INSERTS  \aR  I TE  STATEMENTS  TO  FILE  IOUT 

C 

1 £ 5 CONTINUE 

RETURN 
C 
C 

C FCFMAT  STATEMENTS 
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SUBROUTINE  OUT 


c 

d 0 1 

FORM  AT  </4X 

eo2 

FORMAT (/4X 

€03 

FORM  AT  (53X 

d 04 

FORMAT  </4X 

1 F 1 1 .3  ,F  15 

€05 

FORMAT  (53X 

€06 

FORM  AT (53X 

€09 

FORM AT(4X, 

€ 1 0 

FORMAT (4  X, 

1 F15.0.F10 

€1  1 

FORMAT ( 3 IX 

€12 

FORMAT (3 IX 

eoo 

FORM  AT  (1  HI 

FOR M AT (/4X, 1 3, I 10 ,F13.3,I8,F12.0,3X,5A4,F14.3,F15.0,F10.3,F11.1) 
4,F14.3,F15.C,F10.3,Fll.l) 

FORM  AT  </4X,  13, 1 10  ,F13.3  , I 8 . F I 2 . 0 , 3X , 5HFLC OR , I 3, 12H  COMPARTMENT, 13, 
0 , F 1 0 ,3  ,F11 .1  ) 

5HFLC0R  • 13 , 12H  COMPARTMENT,  13  ,F1 1 • 3 ,F1  5 . 0 , F 1 0 .3  , FI  1 .1  ) 
17HCUTSIDE  DIRECTION,  I3,F14.3,F15.0„F10.3,F11.1> 

3,  F 1 0.  3 »F  1 1 • 0 ) 

FORMAT  (4X,I3,F10.3,F11 . 0 , 3X  , 5HFLOOR  , I 3 , 1 2 H COMPARTMENT,  I3.FH.3, 

FI  1 . 1 ) 

IFLOOR,  13 , 12H  COMPARTMENT  ,I3,F11.3,F15.0,F10.3,FI1.1) 
HOUTSIDE  DIRECTION  , I 3 ,F  1 4 . 3 ,F  1 5 . 0 ,FI  0 . 3 , F 1 1 . 1 ) 

FORMAT (1  HI  , 20X , 18A4  »/S4X  » 8H  A D J LSTED/35X , 4 HTEMP  » 7X , 5HF IXED,28X« 

1 I 2H DIFFERENT I AL, 5X.4HFL0W ,8X , 4HFL0W/4X, 5HFL00R.2X, 1 1 HCQMPARTMENT 

2 » 2X  » 8 F PRESSURE , 2X , 7H PROFILE  » 5 X » 4HFL0W , 3X , 1 fcH CONNECT  ION  TO  , 

E,4X , 1 1HC0EFFICIENT ,2X ,8H  AREA  ,5X,4HFL0W  /) 


3 12  X * 

£01 

FORM, 

£02 

FORM 

£03 

FORM 

£04 

FORM 

1 F 1 1 

£05 

FORM, 

£06 

FORM 

£07 

FORM, 

£C8 

FORM 

l 23F 

2 2ex 

3 5 X , 

4 ,2 X , 

£09 

FORM 

£ 1 0 

FORM 

1 F 15 

811 

FORM 

£12 

FORM 

£13 

FORM 

814 

FORM 

9 0 0 

FORM 

1//5X 

2/5X, 

SOI 

FORM 

I .F15.1.F10.4 


FLOW 


13.1 , I8,F12.0,3X ,5 HFLOOR ,13, 12H  COMPARTMENT, 13, 
Fll.l  ) 

00  R , 1 3 • 1 2H  COMPARTMENT,  1 3 , F 1 1 . 1 ,F1  5 . 1 , F 1 0 .4  , FI  1 . 1 ) 
UTS  IDE  DIRECTION. I3,F14.1,F15.1,F10.4,F11.1) 

1.4H  NET) 

£A  4//20X , 2 OH TEMPERATURE  PROFILE  ,13/  20X, 
COEFFICIENT  , F 1 0 .0//7  2X , 8 HADJ USTED/24X , 5HFI XED , 


€h  AREA  , 5X , 4HFL0W  /) 


TO, 12X  .8HPRESSURE.4X ,1 1 H COEFFICIENT 


COMPARTMENT, 13. FI  1.1 , 


1 //5X , 'FLOW 

2 /5X  ,* PRESSURE 
END 


5HFLC0R , 13 , 12H  COMPARTMENT ,I3.F11.1,F15.1,F10.4,F11.1) 
1 7H0UT  S IDE  DIRECTION  , 1 3 ,F  1 4. 1 ,F1 5 . 1 ,Fl  0 .4  , FI  1 .1  ) 

F8.1  ,4  H NET  ) 

20X  , 18A4) 

iX.'TFE  FOLLOWING  UNITS  ARE  LSED  FOR  OUTPUT* 

IN  LITERS  PER  SECOND  AT  21  DEG  C AND  1 ATM* 

IRE  IN  PASCALS  */EX, ‘AREA  IN  METERS  SQUARED') 

5X , ' THE  FOLLOWING  UNITS  ARE  USED  FOR  OUTPUT* 

IN  CFM  AT  70  DEG  F AND  1 ATM* 


IN 


INCHS  H20 ' / 5X , • ARE  A IN  FEET  SQUARED*) 


SHDG , P 


SUERCUTINE  UNITS. L.l 
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SUBROUTINE  UNITS 


5NBS*PLI ES .SHOW  A. UNITS 
SUERCUTINE  UNITS 
C 

C THIS  ROUTINE  CONVERTS  VARIABLES  H ,FF , Al * AO. CS  TO  SI  UNITS 

C 

c 

PARAMETER  (MM=  1 40 . MS= 8 , MC=9 , MFC=2 , MTP=2 , MFL=25 , M8=50 ) 

COMMON  NT.  P( MM ) ,C( MM.MC) ,NC (MM ) , JC( MM, MC) . ITS( MS ) , 

1 FC(MM.MC), PZ ( MM, MC) »PG( MM.MPO) . CO ( MM , MPQ ) , F ( MM ) , PFO ( MFL , MPO ) , 

2 FF  ( MM  ) , FO ( MM , MPO ) ,CS ( N S ) , PS ( MFL ) ,NS1  ( MS)  ,NS2  < MS)  , 

3 FSS(MS)  ,N ,NS  ,NPO  , ICONV, E, I BUG, AI (MM, MC) ,A0( MM, MPO) ,TITSH(MS,5) 

4 NH,M  MFL)  . IFLOCR  ( MM)  , T ( MTP  , MFL  ) ,NFS1  (MS)  ,NFS2(  MS  ) , I T ( M B ) ,N  TP 

5 ,NCC (MM ), JOC (MM, MPO) .TOUT 
COLELE  PRECISION  P.PO.PS 

C IMEFS  ION  B( 5 ) 

CATA  B/0.3048 .248 .8,0.4719,0.02992,0.0929/ 

CO  10  1=1, NH 

10  F(  I)=H(I  >*B(1  ) 

CO  2 C 1=1, NT 
FF ( I } = FF ( I )*B( 2) 

DO  1C  J=  1 , MC 

AI  ( I ,J  ) = A I ( I , J )*B  ( E) 

16  CONTINUE 

DO  18  1 , MPO 

AO (I ,J  ) = AO ( I, J)*B( E) 

18  CONTINUE 

20  CONT  1NUE 

DO  22  IS=1 ,NS 
22  CS(IS)=CS(IS)*E(4) 

RETURN 

END 


£FDG,P 
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a r.  A n o r>  r>  r>  a a a n r»  A r>  a A a 


SUBROUTINE  SIMEG 


$NBS*PLIB$  .SHOW  A.SIMEQ 
SUBROUTINE  SIMEG 
C 

CFCLESKY'S  METHOD  CF  SOLUTION  OF 
SIMULTANEOUS  LINEAR  ALGEBRIC  EQUATIONS 

PARAMETER  (MM= 1 4 0 ,MS=  e . MC=9 . MFC=2 . MTP=2 . MFL=25  * MB=50 ) 
PARAMETER  (MBF=ME+1) 

COUELE  PRECISION  A.X 
COMMON  /MAT/  A(MB,MBPJ.X(MB) ,N 
NP 1 = N + 1 
2ERO= 1 .OE-35 
K = 0 


SEE  IF  A ( I • 1 ) IS  ZERO 
IF  SC  ADC  ANOTHER  ROW  TO  ROW  I 
IF  (AES  (A  (1 , 1 ) ) .GT.  ZERO ) GO  TO  40 
DO  31  I=1,N 

IF<A(I.l)  a NE  • 0 o ) GO  TO  22 

31  CONTINUE 

12  WRITE(e,804)K 
ST  C F 

32  CC  32  J= 1 , NP 1 

33  A ( 1 , j )=A  (1  . J)+A  ( I , J ) 

CALCULATE  UPPER  AND  LOWER 
TRIANGULAR  MATRICES  OVER  ORIG 
MATR  IX  A 
40  AA  = A (1.1) 

CO  2 J=2,NP1 

2 A( 1 *U)=A(1»J) /A A 
DO  1C  1=2.  N 

K=0 

STORE  A(I.l)  ...  A(I.I)  IN  X ARRAY 
IN  CASE  NEW  A(I. I)  IS  ZERO 
ROW  I CAN  BE  RECALCULATED 

4 CO  5 J = 1 . 1 

5 X(  J ) = A ( I ,J) 

K=K+  1 

CO  10  J=2. NP1 
IF  ( J .GT.  I)  GO  TO  8 
JM 1 = 1 

A A = C . 

DO  3 IR= 1. JM1 

3 AA  = AA+A(  I.  IR)*A(IR. J) 

A ( I ) =A ( I . J ) — AA 

CHECK  IF  A(I, I)  IS  ZERO 

IF  SC  MULTIPLY  OLD  ROW  I BY  2. 

I F ( 1 .NE.  J ) GO  TO  10 

IF ( AE  S ( A ( I » I ) ) .GT.  ZERO) GO  TO  10 

CO  6 JJ=1,I 

6 A(  I,  JJ  ) = X(JJ) 

DO  7 J J=  1 » NP1 
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SUBROUTINE  SIMEQ 


7 A(  I. w)*2.*A( I. J ) 

IF  ( K . CT  . 3)  GO  TO  12 
GO  TC  4 

8 IM  1=1-1 
AA  = C . 

CO  9 IR=  1 « I Ml 

9 AA=AA+A( I,IR)*A(IR. J» 

A< I,  J )=< A< I, J )-AA ) /A (1*1) 

10  CONTINUE 

C ENC  CF  CALCULATION  OF  TRIANGULAR  MATRICES 

C 

C 

EACKkARD  SUBSTITUTION 

X(N)=A<N iNP 1 ) 

00  2 C 1 1 =2  * N 
AA  =0  . 

I=NF  1-  II 
IP  1 = 1+1 
00  IE  J=  IP  1 »N 
15  AA=AA+A ( 1 1 J)*X( J) 

20  X( I ) = A ( I iNPl ) -A A 
C 

804  FORMAT (////// 10X, 1CHPR0GRAM  FAILURE  ,13/////) 
END 


ihDGiP 


FUNCTION  F LO W • L ♦ 1 
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FUNCTIGN  FLOW 


aNES*PLI ES  .SHOW  A.FLCW 

FUNCTION  FLOW (PI.PJ.PZ ,C » 

COUELE  PRECISION  PI.FJ 

THIS  FUNCTION  CALCULATES  FLOWS  EETWEEN  TWO  PO INTS 
C 

IF  C C .LT.  0«001 J G 0 TO  10 

CP=PU-FI+PZ 

SI CN= 1 .0 

IF  ( D F .LT.  • 0 ) S I GN=—  1 . 

FLCW=SIGN*C*SQKT( SIGN*OP) 

RETURN 

10  FLOW  = 0.0 

RETUFN 
END 


a HOG  »P 


FUNCTION  PFLOW.L.l 


r\  r\ 


FUNCTION  PFLOW 


SNBS4PLI  ES • SFC  W A. PFLOW 

FUNCTION  PFLOW  ( I*  P I ) 


C THIS  FUNCTION  CALCULATES  NET  FLOWS  INTO  POINT  I 

C 

PARAMETER  ( MM=140 , M S=e , MC=9 . MFC=2 . MTP=2 * MFL= 25 , MB=50 ) 

COMMON  / CORR/C 1 (M M , MC ) , C2 (MM , MC ) , CO 1 < MM , M PO ) , C02<MM,MP0 ) 

COMMON  NT,  P< MM ) ,C(MM,MC)  ,NC  (MM)  , JC ( MM.MC)  , ITS< MS) , 

1 FC<  MM  ,MC)  ,PZ< MM, MC) , PC ( MM, MPC ) , CO < MM.MPO ) ,F< MM ) , PFO< MFL * MPO ) * 

2 FF ( MM  )»FO(MM»M PO ) * C S ( M S > , PS ( MFL ) ,NSi  < MSI  ,NS2(MS)  * 

3 F SS (MS ) »N»NS , NFO , I CON V , E , I EUG , AI ( MM, MC ) , AO( MM , MPO ) , TI T SH (MS ,5) ♦ 

4 Nh, M MFL) , IFLOOR ( MM) , T( MTP.MFL) ,NFS1 < MS) ,NFS2(MS ),IT(MB),NTP 

5 , NCC< MM  ), JOC( MM, MPC)  ,TOUT 
COUELE  PRECISION  P,PO,PS,PI 
NN=NC ( I) 

SUM=C. 

IF  ( N M .EG.  0 ) GO  TO  3 
CO  1 J J=  1 , NN 
J=  JC  ( I , J J) 

CC  = C1(  1. JJ) 

IF ( F I .LT.  P( J) )CC=C2<  I.JJ) 

PZZ=FZ  (I ,JJ) 

IF  (I  . C-T.  N)PZZ=0. 

FC ( I , J J )=FLOW(PI,P( J ) .PZZ.CC) 

1 SUM=SUM+FC( I , J J ) 

2 NNC=NCC ( I ) 

IF ( NNO  .EQ.  0 ) GO  TO  4 
CO  2 K=l.NNO 
CC  = CO  1(1  , K ) 

IF  ( F I .LT.  PO  ( I «K  ) )CC=CG2  ( I , K ) 

FO ( I ,K  )=FLOW( PI ,PO(I,K) ,0,CC) 

2 SUM=  SUM  + FO ( I , K ) 

4 FFLOH=SUM+FF( I ) 

IF  ( I .LE.  N ) F ( I )=  SUM+FF  ( I ) 

RETURN 

ENC 

aFDG.P  FUNCTION  SFLOW.L.l 
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ft  n o ►-  <n  ft  ft  ft  ft  ft 


FUNCTION  SFLCW 


SNBS  + PLI E$  .SNOW  A.SFLUW 

FUNCTION  SFLOW ( IS  * PI ) 

< 

THIS  RCUTINE  CALCULATES  NET  FLOW  INTO  A SHAFT  AND 
SI-AFT  PRESSURE  PROFILE 


PA R AN E TER  ( M M = 1 4 0 , NS-=8  * MC=9  * MPC=2  , MTP=2 * MFL=25  , MB=50  ) 

COMMCN  NT,  P( MM  ) ,C( MM  ,MC)  ,NC(MM) , JC< MM, MC> . ITS(MS), 

1 FC(NM,MC),PZ(MM,MC), PO ( MM.MPO ) ,CO( MM ,MPO) ,F( MM) , PFO( MFL, MPO ) * 

2 FF  < NM ) ,FO(MM,MFO)  ,CS(  NS) , PS (MFL) ,NS1  (MS) ,NS2(MS) , 

3 FSS (MS  > ,N,NS  ,NPO.  ICCNV , E , I BUG . AI ( MM , MC) , AO(MM.MPO) , T ITSH(MS ,5 ) , 
A NF , F<  NFL) , IFLCOR ( NM) , T ( MTP ,MFL ) ,NF SI ( MS)  ,NFS2( MS) , IT(MB) ,NTP 

5 ,NCC ( NM ), JGC< NN, MFC) .TOUT 
COUELE  PRECISION  P.PO.PS.PI 
IF ( I £ UG  . GT.  1 )WRITE(6.800)IS 
SUM=C. 

N1 =NS1 ( IS) 

N2  = N S2 ( I S ) 

FS(  I ) = PI 
FUF  = C . 

CSS=CS( IS) 

CO  1C  I=N1.N2 
1 1 = I 4 1 — N 1 

FLG=FFLOW<  I,PS( II  ) ) 

FUP  = FLC4-FUP 

SUM=«UM+FLO 

IF ( I .EQ.  N2)GC  TO  5 

IIP1=I  I+l 

SI GN  = 1 

IF ( FLP  .GT.  0 . ) SI GN  =— i . 

FS ( I IP  1 )=PS( 1 1 )-PZ ( I, 1 ) +SIGN*FUP*FUP/( CSS*CSS ) 

IFCIEUG  .GT.  1 ) WR I TE( 6,801)1,11  ,PS( II) ,FLO, FUP*  SUM 
0 CONTINUE 

FS  S ( IS)  = SUM 

SFLO\»=SUM 

RETURN 

FCRNAT  STATEMENTS 

600  FORMAT(///5X,17HFLCW  - SHAFT  NO  ,15/) 

601  FORM  AT (5X.3HI  =,I3.5X,4HII  =,I3,5X.4HPS  =. 

4-  E15.7.5X,  5HFL0  = , E 1 0 . 4 , 5X  , 5HFUP  = , El  0 .4  , 5X,  5HSUM  =,E10.4/) 

END 

aERKPT  PRINTS 
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FEDERAL  INFORMATION  PROCESSING  STANDARD  SOFTWARE  SUMMARY 


01.  Summary  date 

Yr. 

Mo. 

Day 

04.  Software  date 

Yr. 

Mo. 

Day 

2 

0 

5 

2 

8 

08.  Software  type 

09.  Processing  mode 

10. 

Application  area 

General 

Specific 

Automated  Data 

Computer  Systems 

Management/ 

| | System 

| | Interactive 

n Support/Utility 

Business 

Smoke  control  in 

PH  Computer  Program 

jx]  Batch 

pq  Scientific/Engineering  [~ 

Process  Control 

buildings 

| | Subroutine/Module 

| | Combination 

| | Bibliographic/Textual  Q 

Other 

11.  Submitting  organization  and  address 

12.  Technical  contact(s)  and  phone 

02.  Summary  prepared  by  (Name  and  Phone) 

John  Klote  921-3387 


05.  Software  title 

A Computer  Program  for  Analysis  of  Smoke  Control 
Systems 


06.  Short  title 


03.  Summary  action 

New  Replacement  Deletion 

Q □ □ 

Previous  Internal  Software  ID 


07.  Internal  Software  ID 


Suppression  & Extinguishment  Group 
Center  for  Fire  Research 
National  Bureau  of  Standards 
Washington,  D.C.  2023*+ 


John  Klote 
921-3387 


13.  Narrative 


Pressurized  stairwells  and  pressurized  elevators  can  he  used  as  a means  of  providing  a 
smoke  free  exit  route  during  fire  situations.  This  computer  program  analyzes  systems 
intended  to  pressurize  stairwells  and  elevator  shafts. 


14..  Keywords 

Air  movement;  computer  programs;  egress;  elevator  shafts;  escape  means;  modeling; 
pressurization;  simulating;  smoke  control;  stairwells 


15.  Computer  manuf’r  and  model 

Univac  1100/82 

16.  Computer  operating  system 

17.  Programing  language(s) 

ANSI  FORTRAN 

18.  Number  of  source  program  state- 
ments 

1756 

19.  Computer  memory  requirements 

20.  Tape  drives 

21.  Disk/Drum  units 

22.  Terminals 

23.  Other  operational  requirements 

24.  Software  availability 

25.  Documentation  availability 

Available  Limited 

G □ 


26.  FOR  SUBMITTING  ORGANIZATION  USE 


In-house  only 

□ 


Available 

□ 


Inadequate 

□ 


In-house  only 

□ 


Standard  Form  185 

1974  July 

U.S.  Dept,  of  Commerce — NBS 
(FIPS.  Pub.  30) 


185-101 


INSTRUCTIONS 


01.  Summary  Date.  Enter  date  summary  prepared.  Use  Year,  Month,  Day  format:  YYMMDD. 

02.  Summary  Prepared  By.  Enter  name  and  phone  number  (including  area  code)  of  individual  who  prepared  this  summary. 

03.  Summary  Action.  Mark  the  appropriate  box  for  new  summary,  replacement  summary  or  deletion  of  summary.  If  this  software  summary  is  a replace- 
ment, enter  under  “Previous  Internal  Software  ID”  the  internal  software  identification  as  reported  in  item  07  of  the  original  summary,  and  enter  the 
new  internal  software  identification  in  item  07  of  this  form;  complete  all  other  items  as  for  a new  summary.  If  a software  summary  is  to  be  deleted, 
enter  under  “Previous  Internal  Software  ID"  the  internal  software  identification  as  reported  in  item  07  of  the  original  summary;  complete  only  items 
01,  02,  03  and  11  on  this  form. 

04.  Software  Date.  Enter  date  software  was  completed  or  last  updated.  Use  Year,  Month,  Day  format:  YYMMDD. 

05.  Software  Title.  Make  title  as  descriptive  as  possible. 

06.  Short  Title.  (Optional)  Enter  commonly  used  abbreviation  or  acronym  which  identifies  the  software. 

07.  Internal  Software  ID.  Enter  a unique  identification  number  or  code. 

08.  Software  Type.  Mark  the  appropriate  box  for  an  Automated  Data  System  (set  of  computer  programs),  Computer  Program,  or  Subroutine/Module, 

whichever  best  describes  the  software. 

09.  Processing  Mode.  Mark  the  appropriate  box  for  an  Interactive,  Batch,  or  Combination  mode,  whichever  best  describes  the  software. 

10.  Application  Area. 

General:  Mark  the  appropriate  box  which  best  describes  the  general  area  of  application  from  among: 

Computer  Systems  Support/Utility  Process  Control 

Management/Business  Bibliographic/Textual 

Scientific/Engineering  Other 

Specific:  Specify  the  sub-area  of  application;  e.g.:  “COBOL  optimizer”  if  the  general  area  is  “Computer  Systems  Support/Utility”;  “Payroll”  if  the 
general  area  is  "Management/Business”;  etc.  Elaborate  here  if  the  general  area  is  “Other." 

11.  Submitting  Organization  and  Address.  Identify  the  organization  responsible  for  the  software  as  completely  as  possible,  to  the  Branch  or  Division  level, 
but  including  Agency,  Department  (Bureau/Administration),  Service,  Corporation,  Commission,  or  Council.  Fill  in  complete  mailing  address,  including 
mail  code,  street  address,  city,  state,  and  ZIP  code. 

12.  Technical  Contact(s)  and  Phone:  Enter  person(s)  or  office(s)  to  be  contacted  for  technical  information  on  subject  matter  and/or  operational  aspects 
of  software.  Include  telephone  area  code.  Provide  organization  name  and  mailing  address,  if  different  from  that  in  item  11. 

13.  Narrative.  Describe  concisely  the  problem  addressed  and  methods  of  solution.  Include  significant  factors  such  as  special  operating  system  modifi- 
cations, security  concerns,  relationships  to  other  software,  input  and  output  media,  virtual  memory  requirements,  and  unique  hardware  features. 
Cite  references,  if  appropriate. 

14.  Keywords.  List  significant  words  or  phrases  which  reflect  the  functions,  applications  and  features  of  the  software.  Separate  entries  with  semicolons. 

15.  Computer  Manufacturer  and  Model.  Identify  mainframe  computer(s)  on  which  software  is  operational. 

16.  Computer  Operating  System.  Enter  name,  number,  and  release  under  which  software  is  operating.  Identify  enhancements  in  the  Narrative  (item  13). 

17.  Programing  Language(s).  Identify  the  language(s)  in  which  the  software  is  written,  including  version;  e.g.,  ANSI  COBOL,  FORTRAN  V,  SIMSCRIPT  11.5, 
SLEUTH  II. 

18.  Number  of  Source  Program  Statements.  Include  statements  in  this  software,  separate  macros,  called  subroutines,  etc. 

19.  Computer  Memory  Requirements.  Enter  minimum  internal  memory  necessary  to  execute  software,  exclusive  of  memory  required  for  the  operating 
system.  Specify  words,  bytes,  characters,  etc.,  and  number  of  bits  per  unit.  Identify  virtual  memory  requirements  in  the  Narrative  (item  13). 

20.  Tape  Drives.  Identify  number  needed  to  operate  software.  Specify,  if  critical,  manufacturer,  model,  tracks,  recording  density,  etc. 

21.  Disk/Drum  Units.  Identify  number  and  size  (in  same  units  as  "Memory” — item  19)  needed  to  operate  software.  Specify,  if  critical,  manufacturer, 
model,  etc. 

22.  Terminals.  Identify  number  of  terminals  required.  Specify,  if  critical,  type,  speed,  character  set,  screen/line  size,  etc. 

23  Other  Operational  Requirements.  Identify  peripheral  devices,  support  software,  or  related  equipment  not  indicated  above,  e.g.,  optical  character  devices, 
facsimile,  computer-output  microfilm,  graphic  plotters. 

24.  Software  Availability.  Mark  the  appropriate  box  which  best  describes  the  software  availability  from  among:  Available  to  the  Public,  Limited  Availability 
(e.g.:  for  government  use  only),  and  For-ln-house  Use  Only.  If  the  software  is  “Available",  include  a mail  or  phone  contact  point,  as  well  as  the  price 
and  form  in  which  the  software  is  available,  if  possible. 

25.  Documentation  Availability.  Mark  the  appropriate  box  which  best  describes  the  documentation  availability  from  among:  Available  to  the  Public,  In- 
adequate for  Distribution,  and  For  In-house  Use  Only.  If  documentation  is  “Available",  include  a mail  or  phone  contact  point,  as  well  as  the  price  and 
form  in  which  the  documentation  is  available,  if  possible.  If  documentation  is  presently  "Inadequate",  show  the  expected  availability  date. 

26.  For  Submitting  Organization  Use.  This  area  is  provided  for  the  use  of  the  organization  submitting  this  summary.  It  may  contain  any  information 
deemed  useful  for  internal  operation. 
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